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ABSTRACT 

In this work we analyze the spectroscopic properties of a large number of Hn regions, ~2600, located in 38 galaxies. The sample 
of galaxies has been assembled from the face-on spirals in the PINGS survey and a sample described in Marmol-Queralto (2011, 
henceforth Paper I). All the galaxies were observed using Integral Field Spectroscopy with a similar setup, covering their optical 
extension up to ~2.4 effective radii within a wavelength range from ~3700 to ~6900A. 

We develop a new automatic procedure to detect Hn regions, based on the contrast of the Ho- intensity maps extracted from the 
datacubes. Once detected, the procedure provides us with the integrated spectra of each individual segmented region. In total, we 
derive good quality spectroscopic information for ~2600 independent Hn regions/complexes. This is by far the largest nearby 2- 
dimensional spectroscopic survey presented on this kind of regions up-to-date. Even more, our selection criteria and dataset guarantee 
that we cover the regions in an unbiased way, regarding the spatial sampling. 

A well-tested automatic decoupling procedure has been applied to remove the underlying stellar population, deriving the main proper- 
ties (intensity, dispersion and velocity) of the strongest emission lines in the considered wavelength range (covering from [O n] Xil21 
to [Sn] /16731). A final catalogue of the spectroscopic properties of these regions has been created for each galaxy. Additional in- 
formation regarding the morphology, spiral structure, gas kinematics, and surface brightness of the underlying stellar population has 
been included in each catalogue. 

In the current study we focused on the understanding of the average properties of the H n regions and their radial distributions. We 
found that, statistically, there is a significant change of the ionization conditions across the optical extent of the galaxies. The fraction 
of Hn regions located in the intermediate ionization range in a classical BPT diagram is larger for the central regions (r < 0.25r e ), 
than in the outer ones. This is somehow expected, if the origin of this shift is due to the contamination of non-starforming ionization 
sources (e.g., AGN, Shocks, post-AGB stars, etc.), that occur more frequently in the center of the galaxies. 

We find that the gas-phase oxygen abundance and the Ho- equivalent width present negative and positive gradient, respectively. The 
distribution of slopes is statistically compatible with a random Gaussian distribution around the mean value, if the radial distances 
are measured in units of the respective effective radius. No difference in the slope is found for galaxies of different morphologies: 
barred/non-barred, grand-design/flocculent. Therefore, the effective radius is a universal scale length for gradients in the evolution of 
galaxies. 

Other properties have a larger variance across each object, and galaxy by galaxy (like the electron density), without a clear charac- 
teristic value, or they are well described by the average value either galaxy by galaxy or among the different galaxies (like the dust 
attenuation). 

Key words, techniques: spectroscopic - galaxies: abundances - stars: formation - galaxies: ISM - galaxies: stellar content 



* Based on observations collected at the Centro Astronomico 
Hispano Aleman (CAHA) at Calar Alto, operated jointly by the Max- 
Planck Institut fur Astronomie and the Instituto de Astrofisica de 
Andalucfa (CSIC). 



1. Introduction 

Nebular emission lines from bright-individual H n regions have 
been, historically, the main tool at our disposal for the direct 
measurement of the gas-phase abundance at discrete spatial po- 
sitions in galaxies. A good observational understanding of the 
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Fig. 1. Flow chart of HIIexplorer, the algorithm developed to select, segregate, and extract the spectra of the individual H n regions 
for the analyzed galaxies. 



distribution of element abundances across the surface of nearby 
galaxies is necessary to place constrains on theories of galactic 
chemical evolution. The same information is crucial to derive ac- 
curate star formation histories of and obtain information on the 
stellar nucleosynthesis in normal spiral galaxies. 

Several factors dictate the chemical evolution in a galaxy, 
including the primordial composition, the content and distri- 
bution of molecular and neutral gas, the star formation his- 
tory (SFH), feedback, the transport a nd mixing of gas, the 
initial mass function (IMF), et c (e.g. lLopez-SanchezI 1201 Ot 
Lope z-Sanchez & Estebanl 1201 Oi and references therein). All 
these ingredients contribute through a complex process to the 
evolutionary histories of the stars and the galaxies in general. 
Accurate measurements of the present chemical abundance con- 
strain the different possible evolutionary scenarios, and thus the 
importance of determining the elemental composition in a global 
approach, among different galaxy types. 

Previous spectroscopic studies have unveiled some as- 
pects of the complex processes at play between the chem- 
ical abundances of galaxies and their physical properties. 
Although these studies have been successful in determin- 
ing important relationships, scaling laws and systematic pat- 
terns (e.g. luminosity-metallicity, mass- metallicity, and sur- 
face brig htness vs. metallicity relations iLeaueux et al.l 1 1979: 
Skillman 1989 1 IVila-Costas & Edmundsl 119921: IZaritskv etail 
1994t iTremonti et al.l 12004): effective yie ld vs. luminosity and 
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circular velocity relations Garnett 2002; abundance gradients 
and the effective radius of disks iDiaz 119891: systematic differ- 



ences in the g as-phase abundance gradients between n ormal and 
barred spirals IZaritskv et aL| [l994: Martin & Rov 1994; charac- 
teristic vs. integrated abundances Moustakas & Kennicutt 2006; 
etc.), they have been limited by statistics, either in the number of 



observed H n regions or in the coverage of these regions within 
the galaxy surface. 

Hitherto, most studies devoted to the chemical abundance of 
extragalactic nebulae have only been able to measure the first 
two moments of the abundance distribution: the mean metal 
abundances of discs and their radial gradients. Indeed, most 
of the observations targeting nebular emission have been made 
with single-aperture or long-slit spectrographs, resulting in a 
small number of galaxies studied in detail, a small number 
of Hii regions studied per galaxy, and a limited coverage of 
these regions within the galaxy surface. The advent of Multi- 
Object Spectrometers and Integral Field Spectroscopy (IFS) in- 
struments with large fields of view now offers us the opportunity 
to undertake a new generation of emission-line surveys, based 
on samples of scores to hundreds of H n regions and full two- 
dimensional (2D) coverage of the discs of nearby spiral galaxies. 

In the last few years we started a major observational pro- 
gram to understand the statistical properties of H n regions, and 
to unveil the nature of the reported physical relations, using IFS. 
This program was initiated with the PINGS survey (Rosales- 
Ortega et al. 2010), which acquired IFS mosaic data of a num- 
ber of medium size nearby galaxies. In Sanchez et al. (2011) 
and Rosales-Ortega et al. (2011), we studied in detail the ion- 
ized gas and Hn regions of the largest galaxy in the sample 
(NGC628). We then continued the acquisition of IFS data for 
a larger sample of visually classified face-on spiral galaxies 
(iMarmol-Oueralto et al.ll201ll hereafter P aper I), as part of th e 
feasibility studies for the CALIFA survey ([Sanchez e tai1l2012h . 
The spatially resolved properties of a typical galaxy in this sam- 
ple, UGC9837, was presented bv lViironen et al.l ( l2012l) . Face-on 
galaxies are more suitable to study the spatial distribution of the 
properties of Hn regions: (i) they are identified and segregated 
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Table 1. Main properties of the analyzed galaxies 



Name 


RA 


Dec 


Class 


z 


B 


B-V 


mod 


M B 


Vrol 


incl 


PA 


r c 


Type 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10 ) 


(11) 


(12) 


(13) 


(14) 


(15 


2MASXJ1319+53 


13:19:35.9 


+53:30:09.8 


Sb 


0.0248 


15.92 


0.72 


35.3 


-19.3 


51 


28.0 


55 


9.7 


L 


NC 


CGCG 07 1-096 


13:00:33.2 


+10:07:47.8 


Sb(r) 


0.0239 


14.69 


0.79 


35.1 


-20.5 


125 


45.0 


220 


10.7 


I/L 


NC 


CGCG 148-006 


07:44:57.4 


+28:55:39.0 


Sc 


0.0234 


15.16 


0.75 


35.1 


-19.9 


180 


40.5 


250 


10.3 


L 


NC 


CGCG 293-023 


02:30:21.5 


+56:47:29.5 


SBb 


0.0156 


15.51 


0.99 


34.3 


-18.8 


325 


19.6 


100 


6.9 


L 


NC 


CGCG 430-046 


23:00:46.2 


+ 13:37:07.9 


Sc 


0.0243 


14.49 


0.86 


35.1 


-20.7 


168 


42.1 


160 


11.0 


L 


NC 


IC 2204 


07:41:18.1 


+34:13:55.9 


Sab 


0.0155 


15.71 


0.93 


34.2 


-18.5 


49 


24.1 


58 


15.5 


I/E 


AGIn 


MRK 1477 


13:16:14.7 


+41:29:40.1 


SBa(r) 


0.0207 


14.83 


0.47 


34.9 


-20.0 


144 


59.0 


100 


7.9 


L 


NC 


NGC 99 


00:23:59.5 


+ 15:46:13.7 


Sc 


0.0177 


13.71 


0.62 


34.5 


-20.7 


79 


33.7 


40 


19.9 


L 


NC 


NGC 3820 


11:42:04.9 


+ 10:23:03.3 


Sbc 


0.0203 


14.97 


0.95 


34.8 


-19.8 


95 


41.6 


205 


8.5 


L 


NC 


NGC 4109 


12:06:51.1 


+42:59:44.3 


Sa 


0.0235 


14.79 


0.98 


35.1 


-20.7 


132 


44.0 


220 


7.5 


E 


NC 


NGC 7570 


23:16:44.7 


+ 13:28:58.8 


Sa 


0.0157 


13.40 


0.81 


34.2 


-20.8 


81 


65.7 


130 


21.6 


I/E 


NC 


UGC 74 


00:08:44.7 


+04:36:45.1 


Sc 


0.0131 


13.51 


0.85 


33.8 


-20.3 


135 


38.3 


130 


24.6 


L 


(s 


UGC 233 


00:24:42.7 


+ 14:49:28.8 


SBbcD 


0.0176 


14.35 


0.62 


34.4 


-20.1 


161 


4.1 


150 


6.9 


L 


NC 


UGC 463 


00:43:32.4 


+ 14:20:33.2 


SABc 


0.0148 


12.93 


0.85 


34.1 


-21.1 


134 


40.0 


40 


20.6 


L 


(s 


UGC 1081 


01:30:46.6 


+21:26:25.5 


Sc 


0.0104 


13.54 


0.83 


33.3 


-19.8 


82 


22.0 


130 


31.2 


L 


(s 
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01:31:26.6 


+ 14:16:39.0 
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0.0149 


14.56 


0.69 


34.1 


-19.5 


74 


21.6 


60 


19.1 


L 


(s 


UGC 1529 


02:02:31.0 


+ 11:05:35.1 


Sc 


0.0155 


13.67 


1.01 


34.2 


-20.5 


166 


45.4 


132 


17.8 


L 


(s 


UGC 1635 


02:08:27.7 


+06:23:41.7 


Sbc 


0.0115 


14.29 


0.96 


33.5 


-19.2 


71 


32.6 


220 


24.0 


L 


NC 
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02:24:24.8 


-02:09:44.5 


SABc 


0.0046 


13.41 


0.83 


31.5 


-18.1 
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49.4 


25 


67.5 


L 


(s 


UGC 3091 


04:33:56.1 


+01:06:49.5 


SABc 


0.0184 


14.67 


0.88 


34.5 


-19.8 
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39.1 
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5.2 


L 


(s 


UGC 3 140 


04:42:54.9 


+00:37:06.9 
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0.0154 


12.84 


0.97 


34.1 


-21.3 
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L 


AC0S 
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13 


Z4.0 
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i 

L 


API' 

Adz 


UGC 41 07 


07:57:01.9 


+49:34:02.5 


Sc 


0.0117 


13.64 


1.06 


33.7 


-20.9 


84 


19.7 


40 


19.2 


L 


(s 


t i f ^ /"' c 1 r\r\ 

UGC 5100 


09:34:38.6 


+05:50:29.9 


SBb 


0.0184 


14.15 


1.02 


34.5 


-20.4 


326 


69.8 


100 


11.1 


I/L 


(s 


UGC 6410 


11 ^ A f\ C f\ 

11:24:05.9 


AC A O f» 

+45:48:39.9 


SABc 


0.0187 


14.32 


0.78 


34.7 


-20.3 


185 


44.2 


100 


16.5 


L 


NC 


UGC 9837 


15:23:51.7 


+58:03:10.6 


SABc 


0.0089 


13.65 


0.57 


33.2 


-19.5 


90 


23.2 


145 


30.8 


L 


(s 


UGC 9965 


15:40:06.7 


+20:40:50.2 


Sc 


0.0151 


14.06 


0.72 


34.2 


-20.1 


137 


22.4 


73 


17.8 


L 


(s 


UGC 1 1318 


18:39:12.2 


+55:38:30.5 


Sbc 


0.0196 


13.27 


0.72 


34.8 


-21.5 


77 


21.0 


72 
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I/L 


(s 




99-S^-^S Q 
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1 nzL 


^ 1 


9 1 A 




J J.o 


/ J 


9A 9 


t/t 


Mr 
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23:08:57.2 


+ 12:02:52.9 


SABc 


0.0163 


13.26 


0.77 


34.3 


-21.0 


112 


27.2 


40 


24.9 


L 


ac; 
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NGC 628 


01:36:41.7 


+ 15:47:01.0 


Sc 


0.00219 


9.35 


0.56 


29.7 


-20.7 


101 


34.9 


100 


122.2 


L 


AQ 


NGC 1058 


02:43:30.0 


+37:20:29.0 


Sc 


0.00173 


11.24 


0.62 


29.9 


-18.5 


115 


19.6 


100 


45.0 


L 


ac: 


NGC 1637 


04:41:28.2 


-02:51:29.0 


Sc 


0.00239 


11.27 


0.64 


30.4 


-18.5 


114 


31.1 


131 


48.7 


L 


ac; 


NGC 3 184 


10:18:16.8 


+41:25:27.0 


SABc 


0.00194 


10.31 


0.58 


30.2 


-20.0 


108 


24.2 


117 


109.3 


L 


AG 


NGC 3310 


10:38:45.8 


+53:30:12.0 


SABb 


0.00331 


11.15 


0.35 


31.2 


-20.1 


158 


31.2 


237 


37.9 


L 


AC! 


NGC 4625 


12:41:52.7 


+41:16:26.0 


SABm 


0.00203 


12.73 


0.57 


30.5 


-17.8 


105 


46.1 


73 


25.5 


L 


AO 


NGC 5474 


14:05:01.6 


+53:39:44.0 


Sc 


0.00098 


11.11 


0.49 


29.2 


-18.0 


40 


50.2 


100 


58.5 


L 


ac; 



Notes: (1): Galaxy name from NED, the NASA/IPAC Extragalactic Database (http://nedwww.ipac . caltech.edu/ I. (2)-(3): Right ascension 
and declination coordinates in J200 Equinox , expressed in units of RA (hh mm ss) and Dec (dd mm ss). (4): Morphologycal types obtained 
from NED and Hyperleda (Paturel et al. 2003, , http : //leda . univ- lyonl . fr I following the Third Reference Catalogue of Bright Galaxies 
(RC3) classification dCorwin et al .11 1994 , http : //vizier .u-strasbg . fr/viz-bin/VizieR?-source=VII/155 I. (5): Redshift values from 
NED. (6): Dust corrected B-band magnitude calculated from Hyperleda. (7): (B-V) colors obtained from Hyperleda expressed in magnitudes. 
(8) Distance modulus extracted from Hyperleda. (9) Mb: Dust corrected B magnitude extracted from Hyperleda. (10): Rotation velocity in 
units of kms -1 from NED. (11): Absolute value of the inclination angle derived from Hyperleda. (12): Position angle values from our own 
analysis as described in the text, units ar e degrees. (13): Ef f ective radius from this work, units are arcsecs. (14): Disks classification (L:Late, 
E: early, I: intermediate) accord ing to the lLauri kainen et aT] d2010l) criteria. (15): Spiral arm classification following the classes proposed by 
Elmeareen & Elmeareen ( 1987), when not avalaible we use the RC3 classification. 



more easily; (ii) their spatial distribution is less prompt to the 
errors induced by inclination effects; (iii) they are less affected 
by dust extinction along the line of sight within the galaxy and 
(iv) it is more easy to identify their association with a particular 
spiral arm. 

In this article we focus on the study of the properties of the 
Hn regions for the face-on spiral galaxies observed so far. In 
Sect.|2]we summarize the main properties of the analyzed galax- 
ies. In Sect. |3]we present the automatic algorithm developed to 
detect, segregate and extract the integrated spectra of the dif- 
ferent Hn regions within a datacube; a comparison between the 



results derived with this method and those provided with other 
published ones is shown in Sect. |4] In Sect. [5] we describe an 
analytical method to define the presence and location of spiral 
arms within a galaxy. The method has been tested and used to 
associate the detected H n regions to different spiral arms and/or 
to the intra-arm region. The main spectroscopic properties of the 
catalogued H n regions and the morphological structure of each 
galaxy are described in Sect. 16. 1 ,T1 and l6. 1 T2l The main results of 
our analysis are included in Sect. [7] where we describe the sta- 
tistical properties of the H n regions (Sect. I7.lt . and their radial 
gradients (Sect. lT2l i. The conclusions of this study are presented 
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in Sect. [8] In the appendix, the publicly accessible catalogues 
of the properties derived for the analyzed HII regions are de- 
scribed in Sect. [A] and an empirical correction to decontaminate 
the [N II] emission on narrow -band Ha images is proposed in 
SectE 

2. Sample of galaxies 

Table[T]lists the sample of galaxies analyzed in the current study, 
including for each object, its name, central coordinates, mor- 
phological classification, redshift, and some additional infor- 
mation that we will describe later. The sample of galaxies has 
been selected from two available datasets: (1) the IFS survey of 
nearby galaxies described in Paper I, which comprises ~85% of 
the galaxies an alyzed here, and (2) galaxi es selected from the 
PINGS survey (Rosales-Ortega et al. 2010), accounting for the 
remaining ~15% of the galaxies. In both cases we selected vi- 
sually classified face-on galaxies, with a relatively unperturbed 
spiral structure. 

The sample is dominated by late-type spirals (31 out of 
38), a ccording to the classification criteria bv lLaurikainen et all 
shown in Table Q] Therefore, we lack the statistics re- 
quired to analyze possible differences in the properties of H n re- 
gions between late/early-type spirals. About 40% of the galaxies 
show evidence of a bar (based on its visual classification listed in 
TableQ]), although only in a few of them this feature clearly dom- 
inates the morphology of the galaxy (e.g. UGC 5 100). Regarding 
the structure of their spiral arms, the sample includes a mix of 
grand-design spiral s (e.g. NGC628), or clearly flocculent ones 
(e.g., UGC 9837 IViironen et all 1201 2h . Although it is by no 
means a statistically well defined sample, we consider that it is 
representative of the average population of spiral galaxies in the 
local Universe. 

Both sub-samples of galaxies were observed using simi- 
lar techniques (IFS), using the same instrument (PMAS in the 
PPAK configuration, iRoth et al]l200l iKelz et al.ll2006l) . cov- 
ering a similar wavelength range (~3700-6900A), with similar 
resolutions and integration times. The data reduction was per- 
formed using the same procedure (R3D, ISanchez|l2006al) . as 
described in Paperl and Rosales-Ortega et al. (2010). The main 
difference is that for the first sample a single pointing strategy us- 
ing a dithering scheme was applied, while for the largest galaxies 
of the PINGS survey, a mosaic comprising different pointings 
was required. This is due to the differences in projected size, 
considering the different redshift range of both samples: for the 
first one corresponds to ~0. 01-0. 025, while for the second one 
corresponds to ~0. 001-0. 003. Therefore, in both cases the cov- 
ered field-of-view (FoV) corresponds to a similar optical size, 
~2 effective radii, in general (the effective radius is classically 
defined as the radius at which one half of the total light of the 
system is emitted). 

The observational setups allow us to cover the optical wave- 
length range, sampling many of the most important emission 
lines for H n regions, from [O n] /13727 to [S n] A673 1 . Details 
on the observing strategy, setups, reduction and main character- 
istics of the dataset are described in Paperl and Rosales-Ortega 
et al. (2010). The final dataset comprises 38 individual dat- 
acubes, one per galaxy, with a final spatial sampling of 1 "/pixel 
for most of the galaxies. The data cubes were created u sing the 
interpolation scheme described in ISanchez et al.l d2012l) . devel- 
oped for the CALIFA survey. Despite of the original fiber size 
(2.7"/fiber), the three pointing dithering scheme allows to in- 
crease slightly the final resolution. The selected offsets, with val- 
ues corresponding to a fraction of the fiber-size, allows to cover 



the gaps between a djacent fibers too. In average natural seeing 
conditions of ~1" (Sanc hez et al]l2007al) . this technique allows 
to pro vide a final spatial reso lution of FWHM~2" for this instru- 
ment (ISanchez et alJl2007bh . The size of the spaxel was selected 
as the largest convient pixel to the sample this resolution ele- 
ment, 2 pixels per FWHM, i.e., l"/pixel. Due to the large size 
of the IFS mosaics of NGC628 and NGC3184, the two largest 
galaxies observed with PINGS, and the fact that they were not 
observed using a dithering scheme for all the pointings , we set 
the spatial sampling to 2"/pixels. In this case, the final resolu- 
tion is larger than the original size of the fiber, due to the inter- 
polation kernel. A rough estimation indicate that the final spatial 
resolution is FWHM-3.5-4". 

On average the physical spatial sampling ranges between a 
few hundreds of parsecs (for the nearest galaxies) to almost 1 
kpc (for the more distant ones). To derive this physical scale 
it is required to adopt a certain distance modulus. Consistently 
with values reported in TableQ] we adopted the distance modulus 
provided by Hyperleda, mod, defined as: 



mod = 5 log(D L ) + 25 



(1) 



and Di is the luminosity distance in Mpc. The derived scale 
can be compared to the physical diameter of a well-known H n 
region in our Galaxy, i.e. the Orion nebula (D~8 pc), or to the 
extent of those which are considered prototypes of extragalac- 
tic giant H n regions, such as 30 Doradus (D~200 pc), NGC604 
D-46 pc ) or NGC5471 (D~ l kpc) as reported bv lOev etal] 
2003) and iGarcia-Benito et al.l (1201 ll) . Thus, given the under- 
sampling in the physical size of the H n regions in our data, we 
cannot use it to derive direct estimates of the optical extension of 
these regions. In other aspects, like the depth, covered extension 
of the galaxy, projected resolution and wavelength coverage, the 
data provided by both samples are very similar. 



3. Extraction of the H n regions 

The segregation of Hn regions and the extraction of the cor- 
responding spectra is performed using a semi-automatic pro- 
cedure, named HIIexplorejQ. The procedure is based on some 
basic assumptions: (a) Hn regions are peaky/isolated structures 
with a strong ionized gas emission, clearly above the continuum 
emission and the average ionized gas emission across the galaxy; 
(b) H n regions have a typical physi cal size of about a hundred 
or a few hundreds of parsecs (e.g . Gonzalez Delgado & Perezl 
l!997HLopez et al.ll201 lUOev et al.ll2003l) . which corresponds to 
a typical projected size at the distance of the galaxies of a few 
arcsec. 

The algorithm requires a set of input parameters: (i) a line 
emission map, with the same world-coordinate system (WCS) 
and resolution as the input datacube (preferentially an Ha emis- 
sion line map); (ii) a flux intensity threshold for the peak emis- 
sion of each H n region; (iii) a maximum distance to the peak 
location for a pixel associated with each H n region; (iv) a rela- 
tive threshold with respect to this peak emission of the minimum 
intensity of each H n region; (v) an absolute threshold of the min- 
imum intensity corresponding to each H n region. All these pa- 
rameters can be derived from either a visual inspection and/or a 
statistical analysis of the Ha emission line map. The algorithm 
produces as an output a segmentation FITS file describing the 



The adopted modulus for each galaxy is included in the final cata- 
logs, described in the ApendixlAl 

1 http : //www . caha . es/ sanchez/HII_explorer/ 
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LRA(arcsec) 



Fig. 2. IFS-based Ha maps derived for three representative galaxies of the sample (color images), together with the detected Hn 
regions shown as black segmented contours. The regions plotted in grey e have been removed from the final catalogue after a 
sanity-check due to low S/N, quality problems in the extracted spectra, and/or affected by vignetting effects. 



Hllexplorer vs. REGION 

200 F ' ' 




200 100 -100 -200 

ARA (arcsec) 



Fig. 3. Comparison between the Hn region catalogue of 
NGC628 obtained by HIIexplorer us ing IFS data (grey con- 
tours, 286 regions) and the (modified) iFathi et a"D (I2007I) cata- 
logue obtained with Region using a traditional Ha narrow-band 
image (blue circles, 180 regions). Only those regions within the 
FoV of the IFS data have been considered. 



pixels associated to each H n region, designated with a running 
index starting with 1 (e.g. the primary H n region ID in this ar- 
ticle), with the zero reserved to areas not associated with any 



Hn region (i.e. regions free of emission or below the absolute 
threshold described above). 

The segregation algorithm is based on a simple iterative pro- 
cedure, summarized in the flow chart shown in Fig.Q] As a first 
step the algorithm looks for the brightest pixel within the emis- 
sion line map. Its location is stored as the peak/central coordinate 
of a new H n region, associated with a certain running index (ID 
number). After this, the adjacent pixels are aggregated to this 
H n region if all of the following criteria are fulfilled: (i) the dis- 
tance to the central/peak pixel is below the selected limit; (ii) the 
flux at the pixel is above the relative threshold with respect to 
the peak emission; (iii) the flux at the pixel is above the absolute 
flux threshold described before. Whenever any of these criteria 
are not fulfilled the aggregation procedure stops, the ID number 
is increased by one, all the aggregated pixels are masked-off, and 
the peak-identification procedure is repeated. The overall proce- 
dure stops whenever no new peak is detected above the selected 
peak-intensity threshold. 

The outcome of the procedure is illustrated in Fig. [2] where 
we show: (i) the input emission line map, in this case the Ha map 
corresponding to a set of galaxies and (ii) the corresponding de- 
rived segmentation map. Despite the simplicity of the described 
algorithm it is clearly seen that (1) it is able to detect all the 
H n regions that can be identified by eye, and (2) it produces a 
reliable segmentation map. The black segmented areas indicate 
those regions with good quality spectra, while the grey ones in- 
dicate those with poor extracted spectra. The actual procedure 
to detect and reject those ones is described later. In addition, the 
procedure provides with a mask where all the H ii regions are 
flagged out. This mask is important to define the areas where it 
is possible to study the diffuse gas emission. 
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A RA (arcsec) A RA (arcsec) 

Fig. 4. Continuum map of four representative galaxies of the full sample (color images), together with a contour plot of the Ha 
intensity map (both in arbitrary units). The solid lines show the location of the identified spiral arms as a result of the fitting 
procedure. Left side panels show galaxies with clearly defined spiral arms, while right side panels show galaxies with poorly 
identified ones. 



There are other publicly accessible packages for the auto- 
matic sel ection/segrega tion of Hn re gions in the literature (e.g. 
HIIphot. IThilker et al J 12005 Region. iFathi et al.ll2007l) . that in 
principle could be adapted for the main purpose of the cur- 
rent study. However, these packages are strongly focused on the 
analysis of narrow-band images, of much higher spatial resolu- 
tion, where the H n regions are clearly resolved. In some cases 
the procedure requires a detailed knowledge of the observa- 
tional procedure (number of frames co-added, ADU of the CCD, 
etc.). For example, HIIphot uses object recognition techniques to 
make a first guess at the shapes of all sources and then allows for 
departure from such idealized seeds through an iterative grow- 
ing procedure. In ess ence, this algorithm is similar to the one 
used by SExtractor dBertin & Arn outs 1996), for the detection 
and segregation of galaxies in crowed fields. We experimented 



with these packages before developing our own code, but we 
did not get any optimal solution. The main reasons were that (1) 
our data have a much coarser resolution than the one provided 
by narrow-band imaging (even from ground-based telescopes); 
(2) reconstructed IFU map have a strong cross-correlated noise 
among nearby interpolated pixels and (3) none of the preced- 
ing codes provide a final segmentation map usable to extract the 
integrated spectra of the H n regions from the datacubes in a con- 
venient way. 

We experimented with the use of HIIexplorer on a Ha 
narrow-ba nd image provided by the SINGS legacy survey for 
NGC 628 (iKennicutt et alJl2003h . A visual inspect i on of the se- 
lected regions with those shown by IThilker et al.l (|2000), indi- 
cates that although we detect similar regions, HIIexplorer tends 
to define regions of mostly equal size. This is expected, since 
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Table 2. Detected HII regions within our sample. Table 3. Derived parameters for the modelled spiral arms. 
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3.31 


0.36 


NGC 4625 





1 


23 
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9.7 


NGC 3310 
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46 


3.93 


4.25 


2.53 


0.33 


NGC 5474 





2 


60 


1.0 


7.6 


NGC 4625 


66 


49 


3.10 


3.11 


1.83 


0.25 














NGC 5474 


122 


121 


95 


2.74 


2.91 


1.13 


0-14Notes: 


(1) Galaxy name 


used along this article; (2) quality flag 



analysis of the spiral arms: 1 = well defined arms, = arms not well- 
defined; (3) number of arms detected with our modelling; (4) A, B, 
C parameter of the spiral m odel, described in Sect. [5] equation [5] by 
Ringermacher & Mead (2009). A is in units of arcsec, B is a dimension- 
less parameter, and C is in units of radians -1 . 



The adjacent continuum for each pixel was derived by averag- 
ing the flux intensity within two similar bands red- and blue- 
shifted 100 A from the center of the initial one. This contin- 
uum intensity is then subtracted from the Ha intensity to derive 
a continuum-subtracted emission line map. The central wave- 
length of all these bands has been shifted to the observed frame 
taking into account the redshift of the object. The separations be- 
tween the filters and the filter widths are large enough to avoid 
any possible error in the derivation of the Ho- intensity map due 
to kinematic shifts. 

However, this Ha intensity map is contaminated with the ad- 
jacent [Nil] emission lines, and it is not corrected for the emis- 
sion of the underlying stellar population (see Appendix[B] for a 
discussion on the topic). A cleaner Ha emission map can be 
derived using emission-line/stellar population decoupling pro- 
cedures (e.g.. [Rosales-Orteg a et al J 120101 : ISanchez et al.l 1201 lb 



Notes: (1) Name of the galaxy used along this article; (2) number of de- 
tected HII regions; (3) number of HII regions with good quality spectra, 
as described in the text; (4) number of HII regions with H/3 emission 
line detected at 3tr significance; (5) median value of the Ho- intensity 
derived from the narrow-band images in 10~ 16 erg s _1 cnT 2 ; (6) stan- 
dard deviation of the Ho- intensity in the same; (7) median value of the 
estimated radii of the HII regions in units of 100 pc (we need to note 
that the size derived by HIIexplorer is an ill-defined parameter for our 
dataset) ; (8) standard deviation of the previous estimated radii, in the 
same units. 



for the spatial resolution the maximum size allowed for each re- 
gion is reached before that imposed by the ratio of local to peak 
intensity. Our code was never meant to provide a particularly 
reliable measure of the projected size, as in our data this param- 
eter is ill-defined. Specifically, Hn regions can be significantly 
smaller than the resolution element size. In Sect. [4] we present a 
quantitative comparison with methods available in the literature. 

Once tested the procedure, we applied it to our IFS data. 
First, we create a Ha intensity map for each object by co-adding 
the flux intensity within a square-shaped simulated filter cen- 
tered at the wavelength of Ha (6563A), with a width of 60 A. 
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iMarmol-Oueralto et alJfeOl it ISanchez et al.ll2012l) . This allow 
us to recover much fainter emission line regions on top of under- 
lying strong absorption features. However, in most cases these 
emission line regions do not correspond with classical Hn re- 
gions and are asso ciated with other ionization processes (e.g. 
iKehrig et alJl20ll , and references there in). On the other hand, 
the adopted procedure resembles as much as possible the classi- 
cal procedure used to detect H n regions, which will allow us to 
make a better comparison with previous results. 

As the main goal of the current study is to extract the spec- 
troscopic properties of the H n regions, we applied HIIexplorer 
adopting the following input parameters for all the galaxies: (i) 
a minimum flux density for the peak intensity of an H n region 
of 2 10~ 17 erg s _1 crrT 2 arcsec -1 ; (ii) a minimum relative flux 
to the peak intensity for associated pixels corresponding to the 
same Hn region of 10% and (iii) a maximum distance to the 
location of the peak of 3.5" (7" for NGC 628 and NGC 3184). 
The maximum distance was selected using and iterative process, 
maximizing the number of detections when compared with a vi- 
sual inspection, and do not allow to segment clearly single H n 
regions. Then, we extracted a single spectrum for each region 
by co-adding all the spectra in the original cubes with the same 
identification index (ID) in the the derived segmentation map. 
The final spec tra are stored i n the so-called row-stacked spectra 
(RSS) format Sanc hez et alJ 12004). which comprises a 2D spec- 
tral image, with the spectrum corresponding to each H n region 
order by rows (each one corresponding to the considered ID), 
and an associated position table which records the barycenter of 
each H n region (based on the Ha intensity). This format allows 
us to visualize individual spectra of H n regions a nd their spatial 
distri bution using standard techniques (e.g., E3D. ISanchez et al] 
l200l . The ID is a unique identification index that will be used 
to identify the H n regions hereafter (including the tables, figures 
and on-line material). 

Table [2] summarizes the results from this analysis. It shows, 
for each galaxy, the number of H n regions detected, the number 
of regions with good quality extracted spectra, the median Ha 
flux intensity as directly measured from the IFS-based narrow- 
band image, and its standard deviation. Due to the different red- 
shift range, there is a wide variance in the median flux inten- 
sity of the Ha emission, much larger than the absolute luminos- 
ity as we will describe later. Since we are selecting H n regions 
which physical size is slightly smaller than our typical resolu- 
tion element, this implies that we are actually not aggregating a 
large number of real H n regions in each detected complex. Our 
comparisons with higher resolution narrow-band images con- 
firms this suspicion. Finally, we have included in the table the 
median radius of the regions, defined as R= y/A/n, where A is 
the area within a region (Rosales-Ortega et al. 201 1). We should 
state clearly here that the physical scale is an ill-defined param- 
eter in our survey, due to two reasons: 1) the coarse spatial sam- 
pling compared to the expected size of H n regions, these can be 
significantly smaller than the resolution element size and 2) the 
adopted procedure to detect and segregate the regions, namely 
the introduction of an angular upper size limit to the continuous 
emission region. Only for the galaxies at lower redshift, the sizes 
of the H n regions are of the order of the expected one, i.e. ~100 
pc. 

In practice, our segregated Hn regions may comprise sev- 
eral classical ones, in particular for the more distant galaxies. 
Detailed simulations on the effect of resolution loss have shown 
us that on average each selected region at z ~0.02 may com- 
prise 1-3 Hn regions from the ones selected from low redshift 
galaxies, z ~0.002 (Mast et al., in prep.). Following Lope z et aD 



(1201 lh . the considered H n regions would have a size of a few to 
several hundreds of parsecs, based on their Ha luminosity, de- 
tailed in Section 17.11 Table Q Thererefore, the results from the 
simulations are expected, due to the typical size of an extragalac- 
tic Hn region and the lose of physical resolution at the higher 
redshifts. i.e. we are selecting Hn-regions and/or Hn aggrega- 
tions (note that throughout this paper we will refer indistinctively 
to these segmented regions as Hn regions). Therefore, our re- 
sults are not useful to analyze additive/integrated properties on 
individual H n regions, like the Ha luminosity function, but are 
perfectly suited for the study of line ratios, chemical abundances 
and ionization conditions. 

In total, we have detected 3107 Hn regions, 2573 of them 
with good spectroscopic information. To our knowledge this is 
by far the largest 2-dimensional, nearby spectroscopic H n region 
survey ever accomplished. 

4. Comparisons with previous selection methods to 
detect H ii regions 

In order to assess quantitatively the degree of segmentation pro- 
vided by HIIexplorer with respect to other traditional H n region 
catalogues generators, we performed a comparison between the 
H n region catalogues of NGC 628, NGC 3 1 84 and NGC 5474 
obtained with HIIexplorer and those reported in the literature. 

For NGC 628 we used the Ha image extracted from the IFS 
data with a resolution of 2 arcsec/pixel and a Fo V of ~ 6 arcmin, 
and th e one produced by the REGION software in Fath i et al.l 
( 2007) (hereafter Fa07), obtained from a narrow-band Ha image 
with a resolution of 0.33 arcsec/pixel and a FoV of ~ 1 1 arcmin. 
Similar comparisons could be performed with an y other pack- 
age c reated to segregate H n regions (e.g. HIIPhot. iThilker et al.l 
2000). We selected this one because Fa07 provided publicly ac- 
cessible catalogues. 

The full Fa07 catalogue has 376 regions of which 299 are 
within the FoV of our IFS data. However, the public Fa07 cata- 
logue reports only the position and the full-width at half maxi- 
mum (FWHM) of each H ii region, not the actual shape obtained 
by the software, and this leads to significant overlaps when the 
REGION catalogue is plotted over the galaxy image. Taking 
this into account and considering the difference in resolution be- 
tween the two Ha images, we created a modified version of the 
Fa07 catalogue in order to make a fair comparison. The modified 
Fa07 catalogue was obtained in an iterative way. First we took 
the first region of the catalogue and calculated its distance from 
the rest of the regions. Those regions for which the distance was 
less or equal to the sum of their radii were considered as a sin- 
gle region. In this case, the involved regions are removed from 
the original catalogue and a new entry is added with coordinates 
and size corresponding to the luminosity weighted mean of the 
merged regions, the process is repeated for the rest of the cat- 
alogue entries in an iterative manner. We obtain 180 regions in 
the modified Fa07 catalogue of NGC 628 . 

Fig.[3]shows the comparison between the modified Fa07 cat- 
alogue and the HIIexplorer segmentation map. The blue circles 
correspond to the modified Fa07 catalogue, while the grey con- 
tours to the 286 segmented regions obtained by HIIexplorer for 
NGC 628 based on the IFS data. We note that, a) HIIexplorer 
detects and segments more regions than Fa07, except for those 
cases in which the difference in spatial resolution (0.33 vs. 2 
arcsec/pixel) prevents further segmentation; b) There is a nearly 
1:1 correspondence of regions detected in Fa07 with respect 
to HIIexplorer, the incompleteness of Fa07 with respect to 
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HIIexplorer is 5%; c) 19% of the regions in HIIexplorer have 
2 or more regions of the modified Fa07 catalogue, which is sim- 
ply due to the difference in resolution. We have checked visually 
the extracted spectra of the additional H n regions detected by 
our algorithm, and inspected the original narrow band image and 
they seem to be real H n regions, clearly distingued from the low 
surface brightness diffuse gas. The performance of HIIexplorer 
compared with REGION is remarkable, considering both that 
the narrow-band Ha image used to generate the Fa07 catalogue 
is deeper than the image extracted from the IFS data, and that 
HIIexplorer runs in a completely automated way. 

Th e Hn regions of NGC3184 and NGC5474 were stud- 
ied bv lBradlev etall (120061) using REGION (hereafter B06). For 
NGC 3184, the catalogue obtained by B06 contains 576 H n re- 
gions of which 209 are within the FoV of our IFS data. Like 
the Fathi et al. (2007) case for NGC628, the B06 catalogues re- 
port only the offset from the galaxy centre and the total area 
of the region, not the actual shape obtained by the software, 
which leads to significant overlaps when the REGION B06 cat- 
alogue is plotted over a RA vs. Dec plane using an effective 
radius derived from the B06 catalogue. Therefore, we applied 
the same methodology for a fair comparison obtaining a mod- 
ified B06 catalogue for this galaxy, imposing a Ho- luminosity 
threshold of log(LH ff ) > 37.96 erg s (the minimum luminosity 
detected by HIIexplorer at this redshift). The level of complete- 
ness is 73%, i.e. regions detected by HIIexplorer with respect 
to B06 (note that in the majority of cases there is a 1:1 corre- 
spondence); in 15 cases 2 or more B06 regions are found within 
1 segmented area by HIIexplorer. However, in 5 cases 2 Hn 
regions by HIIexplorer correspond to 1 region found by B06, 
while 13 regions detected by HIIexplorer are not present in the 
B06 catalogue. 

In the case of NGC 5474, the original B06 catalogue con- 
tains 165 Hn regions, of which 98 are within the FoV of the 
IFS data. For this galaxy, we worked directly with the published 
catalogue without further modifications for a better visual com- 
parison. There was no need to apply a luminosity threshold since 
all the regions were above the minimum luminosity observed by 
the regions segmented by HIIexplorer, log(Le ff ) ~ 36.6 erg s" 1 . 
We note that HIIexplorer detects and segments more regions 
than B06, except for those cases in which the difference in spa- 
tial resolution prevents further segmentation. The level of com- 
pleteness (regions detected by HIIexplorer compared to the B06 
catalogue) is of 90% (including 1 : 1 correspondence and multi- 
ple B06 HII regions within one HIIexplorer segmentation), but 
interestingly 3 1 regions detected by HIIexplorer are not found 
in the B06 catalogue, which is surprising given the that the Ha 
image used to generate the B06 catalogue is deeper than the Ha 
map extracted from the IFS data. 

This exercise shows that HIIexplorer is capable of perform- 
ing an excellent H n region extraction for the resolution of our 
IFS data, and that the generated catalogues are comparable (and 
even more efficient) than those generated in a traditional way 
based on narrow-band Ha imaging. 

5. Analytical characterization of the spiral arms 

A fundamental question regarding the star-forming regions in 
galaxies is whether their distribution and properties depend on 
their association (or not) with a particular spiral arm. Two main 
questions are directly connected with this one: (i) whether there 
are azimuthal variations within the spectroscopic properties of 
the H n regions, which would possibly reflect non-radial differ- 
ences in the galaxies evolution, maybe induced by non secular 



processes, and (ii) whether the properties of the inter-arm H n 
regions are different than those of the H n intra-arms ones, which 
will reflect a possible differential evolution associated with ram 
pressure in the spiral arms. The lack of a sample with a statisti- 
cally large number of H n regions, with homogeneously derived 
spectroscopic properties, and with a good characterization of the 
structure of the spiral arms has not allowed to give a conclusive 
answer to these questions so far. 

In the following we attempt to give a good description of the 
structure (number, shape, radial path) of the spiral arms and to 
define a procedure to associate H n regions with each spiral arm 
and/or classify t hem as inter-arm ones. We ado pt the prescrip- 
tion proposed by Ringermac her & Mead! (Hp09) to describe the 
general shape of the spiral arms. This formalism describes the 
radial path of any spiral arm using the formula: 



log(fitan(6»/2C)) 

This function intrinsically generates a bar in a continuous, 
fixed relationship relative to an arm of arbitrary winding sweep. 
A is simply a scale parameter while B, together with C, deter- 
mines the spiral pitch. Roughly, larger C results in tighter wind- 
ing. Greater B results in larger arm sweep and smaller bar/bulge, 
while smaller B fits larger bar/bulge with a sharper bar/arm junc- 
tion. Thus B controls the "bar/bulge-to-arm" size, while C con- 
trols the tightness much like the Hubble scheme. Special shapes 
such as ring galaxies with inward and outward arms are also de- 
scribed by the analytic continuation of the same formula, which 
is particularly useful to analyze the diversity of spiral structures 
within our sample. 

The previous formula describes the radial path in the phys- 
ical plane of the disk for each spiral arm. To describe the full 
observed spiral structure, for a galaxy with Na arms, it is re- 
quired to project it at the observed plane (taking into account the 
inclination and position angle), and to add Na copies of the con- 
sidered arm, rotated by an angle of 360/Wa degrees with respect 
to the precedent one. 

The optimal parameters that describe the current spiral struc- 
ture (Na, A, B and C) have been derived using an interactive 
fitting algorithm that is based on two simple assumptions: (1) 
the spiral arms trace the location of the stronger Ha emission 
and continuum emission, i.e., the integrated intensity along the 
arm should be maximized; and (2) the H n regions are more fre- 
quently clustered around the spiral arms, i.e., the distance from 
each region to the nearest spiral arm should be minimized. The 
analytical parameters of the spiral arms are then derived based 
on these assumptions and using as inputs (i) a broad band im- 
age of the galaxy. The SDSS g-band image i n most of the case s 
(extracted from the SDSS imaging survey, I York et all 12000. , 
and Paper I), and when not feasible the V-band one, from our 
own observations (Paper I), or literature data (Rosales-Ortega et 
al. 2010); (ii) the spatial distribution of Hn regions derived by 
HIIexplorer; and (iii) a few simple assumptions of the number 
of arms and the scale-length of the possible bar and/or the initial 
ring, based on the visual inspection of the images. In general, we 
selected the spiral structure with the smallest possible number of 
spiral arms that fulfill the criteria. 

We are aware that the formalism adopted here to describe 
the analytical structure of the spiral arms is clearly not the most 
mathematically exact one. However it is useful for the ultimate 
goal of our study, i.e. to determine how many spiral arms are in 
a considered galaxy and if a certain H n region belongs to an arm 
or not. A more analytical description of the spiral arms is clearly 
out of the scope of the current study. 
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Table 4. General properties derived for the 10 brighest HII regions within UGC9837. 



ID 


RA 


Dec 


X bs 


Y bs 


x rcs 


Y rcs 


R 


6 


N fl 


F 


D;n in 


D sp 




Vrot 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


UuL yaj /-UU1 


Tin OAon 




Z4.U 


ty.i 


t).4 


0.0 


0.8 


00.3 


1 


i 
1 


l.J 


Z44.Z 




zooj.U 


jy. ij 


UGC 9837-002 


230.9611 


58.0549 


-29.7 


6.9 


4.6 


-4.5 


6.5 


315.6 


4 


1 


1.9 


123.2 


-3.0 


2715.1 


39.56 


UGC 9837-003 


230.9636 


58.0579 


-12.5 


17.9 


4.5 


0.0 


4.5 


0.0 


4 


1 


6.0 


29.5 


-13.4 


2723.2 


39.42 


UGC 9837-004 


230.9661 


58.0574 


4.6 


16.0 


2.1 


2.9 


3.6 


53.3 


3 


1 


0.9 


66.2 


2.9 


2680.2 


39.35 


UGC 9837-005 


230.9650 


58.0580 


-2.5 


18.1 


3.3 


1.8 


3.8 


29.1 


3 


1 


4.7 


121.9 


4.6 


2694.1 


39.24 


UGC 9837-006 


230.9670 


58.0479 


10.9 


-18.3 


-4.3 


-0.3 


4.3 


184.6 


1 


1 


3.0 


220.5 


5.9 


2602.2 


39.17 


UGC 9837-007 


230.9673 


58.0559 


13.1 


10.8 


0.3 


3.7 


3.8 


85.8 


3 





7.1 


35.0 


-20.5 


2639.4 


39.27 


UGC 9837-008 


230.9625 


58.0540 


-20.1 


4.0 


3.0 


-3.1 


4.4 


313.8 


4 





6.4 


78.4 


11.7 


2702.2 


38.96 


UGC 9837-009 


230.9623 


58.0588 


-21.0 


21.1 


6.0 


-1.1 


6.1 


349.3 


3 


1 


6.3 


220.5 


7.9 


2737.6 


38.99 


UGC 9837-010 


230.9644 


58.0596 


-6.6 


24.0 


4.8 


1.9 


5.1 


21.3 


3 


1 


1.6 


143.2 


3.8 


2712.7 


39.13 



Notes: (1) Unique ID of the HII region; (2) Right ascension of the HII region, in degrees, for the J1200 equinox; (3) Declination of the HII region, 
in degrees, for the J1200 equinox; (4) Relative right ascension from the center of the galaxy, in arcsecs; (5) Relative declination from the center of 
the galaxy, in arcsecs; (6) Deprojected X coordinate from the center, in kpc; (7) Deprojected Y coordinate from the center, in kpc; (8) Deprojected 
distance to the center, in kpc; (9) Deprojected position angle, in degrees; (10) ID of the nearest spiral arm, used to associate a HII region with a 
particular arm; (11) Flag indicating if the HII region is clearly associated to the corresponding spiral arm (1) or not (0); (12) Minimum distance 
to the nearest spiral arm, in arcsec; (13) Spiralcentric distance, i.e. distance along the nearest spiral arm from the center, in arcsec; (14) Angular 
distance to the nearest spiral arm, in degrees; (15) Velocity of the ionized gas derived the fitting to the Ha emission line of the HII region, in km/s; 
(16) Decimal logarithm of the dust corrected absolute luminosity of Ho-, in units of erg s _1 . The full catalogue of Hn regions for this object and 
the remaining ones discussed along this article are listed in Appendix lAl including the errors of the velocity and the Ha luminosity. 



Table [3] lists the results of this analysis, including, for each 
galaxy, its name, the number of derived spiral arms, a flag in- 
dicating the reliability of the results, and the parameters of the 
radial path for each arm, according to the described formula: A, 
B and C. The quality flag is 1 for those galaxies with a clearly 
distinguised spiral structure and a well defined set of parame- 
ters to describe them, and for those galaxies without a well de- 
fined spiral structure (flocculent). Fig.|4]shows four examples of 
the derived spiral structure for the galaxies in our sample, in- 
cluding two cases with well defined spiral structure (NGC628 
and UGC 1081), one without a clear defined spiral structure 
(NGC 1637) and a clear flocculent case (UGC 9837). 



We associated H n regions to the nearest spiral arm by com- 
puting the minimum distance between the centroid of the H n re- 
gion and the radial distribution of the considered arm. The mean 
of these distances is then used as a scale-length to separate be- 
tween H ii regions clearly associated with an arm, and possible 
inter-arm ones. A final flag is included in the corresponding cat- 
alogue table describing the coordinates of the detected H n re- 
gions, indicating the nearest arm and the relative distance with 
respect to the median one. Table |4] illustrates the result of this 
analysis. It shows, for one single galaxy (UGC 9837), the abso- 
lute, relative, polar and deprojected coordinates of the 10 bright- 
est H n regions, together with identification of the nearest spiral 
arm, the Cartesian and angular distance to this arm and the spi- 
ralcentric distance (i.e. the distance to the center along the spiral 
arm). In this context brightness means the peak intensity within 
a certain Hn, as considered by HIIexplorer . In addition, we 
included in this table the systemic velocity and absolute lumi- 
nosity of Ha for each H n, derived on basis of the emission line 
fitting described in Section lo*. 1 . 1 1 Similar parameters are derived 
for all the H n regions in the different galaxies, as indicated in 
Appendix [A] To our knowledge, this is the first attempt to per- 
form an analytical association of H n regions to a particular arm 
and/or to an inter-arm area in a survey mode. 



6. Deriving the main spectroscopic properties of 
the H ii regions 

6.1. Decoupling the emission lines from the underlying 
stellar population. 

To extract the nebular physical information of each individual 
H ii region, the underlying continuum must be decoupled from 
the emission lines for each of the analyzed spectra. Several 
different tools have been developed to model the underlying 
stellar pop ulation, effectively decoupling i t from the emission 
lines (e.g.. [Ca pnellari & Emselle mj l2004t [Cid Fernandes et all 
2005l:IOcvirk et al l2006USarzi et all2006l:ISanchez et al.l2 007a: 
Koleva et al.ll2009l: [MacArthur et al.ll2009t IWalcher et al.ll201 U 
Sanchez et al.ll201 lb . Most of these tools are based on the same 
principles, i.e., they assume that the stellar emission is the result 
of the combination of different (or a single) simple stellar pop- 
ulations (SSP), and/or the result of a particular star-formation 
history, whose emission is redshifted due to a certain systemic 
velocity, broadened and smoothed due a certain velocity disper- 
sion and attenuated due to a certain dust content. 

For the particular case of the H n regions, the main purpose 
of this analysis is to provide a reliable subtraction of the under- 
lying stellar population. For doing so, we performed a simple but 
robust modelin g of the continuum emission. We use the routines 
describ ed in ISanchez et al.l d2011l) and iRosales-Or tega et al. 
(l2010h . which provided us with certain parameters describ- 
ing the physical components of the stellar populations (e.g., 
luminosity-weighted ages, metallicities and dust attenuation, to- 
gether with the systemic velocity and velocity dispersion) and 
a set of parameters describing each of the analyzed emission 
lines (intensity, velocity and velocity dispersion). A simple SSP 
template grid was adopted, consisting of three ages (0.09, 1.00 
and 17.78 Gyr) and two metallicities (0.0004 and 0.03). The 
models were extracted from the SSP templa te library provided 
by the MILES project (Vazd ekis et al.ll2010h . The two consid- 
ered metallicities are the most metal poor and most metal rich 
with the largest coverage range in ages, within the considered 
library. The oldest stellar population was selected to reproduce 
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Fig. 5. Spectra of three typical Hn regions for two different galaxies, NGC628 (left panels), and UGC 1081 (right panels). The 
panels correspond, from top to bottom, to the brightest, average and faintest H n regions with good quality data for each particular 
object. The spectral range shows all the emission lines analyzed in this article . From top to bottom: each panel shows the integrated 
spectra of the considered H n region (solid-black), together with the best model for the underlying stellar population (solid-red), and 
the pure-gas spectra (solid-blue). The two boxes show a zoom of the same plots at the Y\j3 and Ha wavelength range. 



the reddest possible underlying stellar population, mostly due 
to larger metallicities than the one considered in our simplified 
model, although it is clearly older than the accepted cosmologi- 
cal time of the Universe. Our youngest stellar population is the 
2nd youngest in the MILES library with both extreme metallic- 



ities. No appreciable difference was found between using this 
one or the youngest one (~80 Myr). Finally, we selected an av- 
erage stellar population, of ~lGyr, required to reproduce the 
intermediate-to-blue stellar populations, and to produce more re- 
liable corrections of the underlying stellar absorptions (Paper!). 
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Table 5. Emission line ratios with respect to H/3, and H/3 line intensity for the 10 brighest HII regions within UGC9837. 



ID 




t(H/j) 


HP 


LOiUJSOO? 

US 


IOIJ6300 
US 




Ha 

H/3 


Hp 


HeI667S 

m 


LS11J67 17 

Hp 


IS11J67J1 
HP 




on 1 


ZU.OU+ 0. IJ 


i on_i_ o o^ 
t .zU± U.Uj 


4.44± 


0.20 


0.08± 


0.02 


A A /I _i_ 

4.44+ 


U.z / 


0.26± 


0.11 


U.Uo± U.U1 


0.34± 


0.03 


0.25± 


0.03 


UGC9837- 


■002 


36.03± 0.66 


2.46± 0.07 


2.07+ 


0.06 


0.08± 


0.01 


3.31+ 


0.17 


0.37+ 


0.12 


0.05± 0.01 


0.48± 


0.03 


0.33± 


0.03 


UGC9837- 


■003 


29.28± 0.63 


2.97+ 0.13 


1.03± 


0.04 


0.08± 


0.02 


3.19+ 


0.21 


0.55± 


0.16 


0.05± 0.01 


0.65± 


0.04 


0.45± 


0.04 


UGC9837- 


■004 


25.17+ 0.57 


3.14+ 0.13 


1.27± 


0.05 


0.11+ 


0.02 


3.18+ 


0.21 


0.55± 


0.15 


0.02± 0.01 


0.62± 


0.04 


0.44± 


0.04 


UGC9837- 


■005 


13.09± 0.51 


3.04+ 0.26 


0.85± 


0.07 


0.11+ 


0.04 


3.57± 


0.31 


0.69± 


0.19 


0.06± 0.01 


0.67± 


0.08 


0.52± 


0.08 


UGC9837- 


■006 


19.40± 0.54 


3.04+ 0.13 


1.77+ 


0.08 


0.08± 


0.04 


3.04± 


0.21 


0.38± 


0.13 


0.02± 0.01 


0.50± 


0.05 


0.37+ 


0.04 


UGC9837- 


■007 


13.90± 0.54 


2.93± 0.19 


1.48+ 


0.10 


0.10± 


0.03 


3.58± 


0.29 


0.50± 


0.17 


0.07± 0.01 


0.61± 


0.06 


0.45± 


0.06 


UGC9837- 


■008 


12.07± 0.45 


3.53+ 0.27 


0.97± 


0.07 


0.13+ 


0.04 


3.04± 


0.28 


0.59± 


0.19 


0.03± 0.01 


0.71± 


0.07 


0.51± 


0.06 


UGC9837- 


■009 


9.31±0.48 


2.10+ 0.14 


3.52± 


0.23 






3.35± 


0.27 


0.27± 


0.11 


0.04± 0.01 


0.38± 


0.06 


0.29± 


0.06 


UGC9837- 


■010 


13.46+ 0.51 


3.52± 0.23 


1.35+ 


0.09 


0.18+ 


0.05 


3.30± 


0.29 


0.52± 


0.19 


0.04± 0.01 


0.73± 


0.07 


0.52± 


0.06 



Notes: H/? fluxes are in units of 10~ 16 erg s _I cirT 2 ; All line intensities have been derived after subtracting the underlying stellar population, 
but without any further correction. The full catalogue of emission line ratios for the Hn regions analyzed in this object and the remaining ones 
discussed along this article are listed in Appendix lAl 

Table 6. Emission line equivalent with of the strongest emission lines analyzed, for the 10 brighest HII regions within UGC9837. 



ID 










Equivalent Width 












[Oil] 


w 


[OIII] 


[OI] 


Ha 


[Nil] 


Hel 


[SII] 






,13727 




,15007 


,16300 




,16583 


,16678 


,16717 + 6731 


UGC9837- 


■001 


-80.6± 22.1 


-83.8+ 17.6 


-367.1+ 208.1 


-8.0+ 1.2 


-503.6+ 497.3 


-29.6+ 33.4 


-6.3+ 0.3 


-69.5+ 7.1 


UGC9837- 


■002 


-102.7+ 40.3 


-53.9± 6.4 


-113.7+ 19.5 


-6.0+ 0.5 


-269.6+ 94.5 


-30.2+ 13.4 


-3.7+ 0.2 


-69.2+ 5.5 


UGC9837- 


■003 


-66.8± 19.1 


-22.7± 2.1 


-24.1+ 1.2 


-2.2+ 0.2 


-96.5+ 11.9 


-16.5+ 3.2 


-1.4+ 0.1 


-34.0+ 1.3 


UGC9837- 


■004 


-66.8± 23.8 


-23.2± 2.2 


-31.0+ 1.9 


-3.2+ 0.3 


-101.6+ 13.4 


-17.9+3.6 


-0.8+ 0.1 


-34.6+ 1.4 


UGC9837- 


■005 


-42.9± 6.7 


-14.0± 1.4 


-12.2+ 0.6 


-1.9+ 0.3 


-64.1+ 7.8 


-12.3+ 2.3 


-1.1+ 0.2 


-22.1+ 0.8 


UGC9837- 


■006 


-76.2+ 29.8 


-30.0+ 2.5 


-54.0+ 4.3 


-3.1+ 0.7 


-129.2+ 31.3 


-16.3+ 5.5 


-0.7+ 0.3 


-38.1+ 1.8 


UGC9837- 


■007 


-51.2+ 10.7 


-21.1+ 2.4 


-31.7+ 2.3 


-2.6+ 0.4 


-97.7+ 18.6 


-13.5+3.7 


-1.8+ 0.1 


-30.2+ 1.3 


UGC9837- 


■008 


-55.0+ 18.5 


-15.2+ 1.7 


-15.1+0.7 


-2.4+ 0.3 


-59.9+ 5.9 


-11.8+ 2.0 


-0.5+ 0.1 


-25.1+0.8 


UGC9837- 


■009 


-40.2+ 7.3 


-22.9± 2.9 


-80.9+ 11.1 


-2.6+ 1.4 


-109.0+ 12.6 


-8.5+ 1.9 


-1.3+ 0.3 


-21.5+ 1.2 


UGC9837- 


■010 


-53.4± 16.2 


-16.5+ 1.7 


-22.7+ 1.5 


-4.0+ 0.5 


-72.0+ 7.9 


-11.4+ 2.3 


-0.9+ 0.2 


-28.1+ 1.4 



Notes: All the listed equivalent widths are in units of A. The full catalogue of equivalent widths for the H n regions analyzed in this object and the 
remaining ones discussed along this article are listed in Appendix lAl 



This library is clearly insufficient to describe in detail the na- 
ture of all the stellar populations and star formation histories. 
However, it covers the parameter space of possible stellar pop- 
ulations well enough to describe them at to 1st order, providing 
a clean residual information of the ionized gas. Evenmore, with 
a combination of the considered templates it is possible to re- 
construct any of the SSP of the full MILES library within an ac- 
curacy similar to our photometric uncertainty (~10%, Paper I). 
Therefore, to include any other template is redundant for the 
main purpose of this analysis. 

Fig.[5]illustrates the results of the fitting procedure for three 
H ii regions (the brightest, the average and the faintest one, in 
terms of Ha luminosity) extracted from two typical galaxies 
(NGC628 and UGC 1081). The figure shows for each Hn re- 
gion the extracted spectrum (black line), together with the best 
multi-SSP model for the stellar population (red line), and the 
pure nebular emission spectrum (blue line). 

6.1 .1 . Deriving the main properties of the emission lines 

To derive the properties of the stronger emission lines detected 
in the stellar-population subtracted spectra, each line was fit- 
ted with a single Gaussian, coupled with the systemic velocity 
and velocity dispersion of different emission lines when needed 
(e.g., for doublets and triplets). This procedure provide us with 



the intensity, systemic velocity and velocity dispersion for each 
emission line. Note that by subtracting a stellar continuum model 
derived with a set of SSP templates, we are already taking into 
account (and correcting for, to a first order) the contribution of 
underlying absorption, which is particularly important in the Ha 
and H/3 lines. 

As discussed in Paper I, there are different issues that may 
affect the final quality of the individual spectra in the datacubes 
(uncleaned cosmic rays, trace problems, low transmission fibers, 
spectra near to the edge of the FoV and vignetting). These qual- 
ity issues, that affect a reduced number of spaxels, propagate 
along the segregation, extraction and "stellar continuum" clean- 
ing processes, and therefore they may affect the final quality 
of the "pure emission" spectra of the Hn regions. In order to 
minimize the impact of these issues on the final sample of H n 
spectra, we have performed an automatic quality check. Only 
those "pure emission" spectra fulfilling the following criteria are 
flagged as good quality data: 

1. The derived intensity for Ha is above zero or below three 
times that of the brightest H n region based on the narrow- 
band image intensity. The contrary may happen in case of 
problems with the fitting procedure, or problems with one 
or a few spectra of those that were co-added to derive the 
integrated spectrum (like a cosmic rays). 
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Fig. 6. Left-panel: Ha velocity field derived from the analysis of the H n regions for UGC 9965. Central-panel: Best model fitted to 
the velocity map using a simple arctan velocity curve, and fitting the inclination, position angle and maximum rotational velocity. 
Right-panel: Residual of the subtraction of the model to the veloc ity map. Each plotted sy mbol in each panel represent an indi- 
vidual H n region. The circles represent those H n regions below the Kauffm ann et al.l (|2003f) demarcation line in the classical BPT 



diagram (Bald win etail l981). and the squares corresponds to those ones located in the intermediate region between that line and 
the Kewlev et al. I (120011) one, as described later in the text. The size of the symbols are proportional to the Ha intensity 



2. The fraction of spectral pixels with negative values in the 
original spectra, i.e. prior to the subtraction of the underly- 
ing stellar continuum, is lower than 10%. The contrary may 
happen in the outer regions of the galaxies, if there is any 
problem with the sky subtraction. 

3. The fraction of spectral pixels in the "pure emission" spec- 
tra of the Hn region with a value below the median flux lcr 
within the wavelength range between 3900 and 6500A, is at 
maximum three times lower than the median of this fraction 
for all the H n regions in the same galaxy. This criterion is 
used to remove spectra strongly affected by the vignetting 
effect, which affects only ~30% of the data (Paper I). 

4. The derived intensity for Ha is more than five times 
above the background noise {aback), estimated as (Tback = 
0"63OO-65OO FWHMu„ e , where (i) o" 6 300-6500 is the standard 
deviation of the continuum intensity once subtracted the un- 
derlying stellar component for the wavelength range between 
6300 and 6500A(i.e., a continuum adjacent to Ha); and (ii) 
FWHMii„ e is the full width at half maximum derived for the 
emission line, as described before. 

The remaining regions are flagged out, masked, and the cor- 
responding spectra are set to zero. The criteria were based on 
iterative experiments on the data, and visual inspections of hun- 
dreds of spectra, before and after subtracting the underlying 
continuum. Although the fraction of flagged-out/rejected spec- 
tra change from object to object, on average this final clean- 
ing affects ~15% of the Hn regions, as can be seen in Table [2] 
Finally, only those H n regions with measured HB emission line 
detected at >3<x significanc e were considered for further analy- 
sis (e.g. lMarino et al.ll2012l) . although they were not masked out 
and their spectra were not set to zero. This criteria was included 
to consider only those regions with good line diagnostic ratios 
and well defined Balmer ratio, both required in further analysis. 
It further reduces the number of selected H n regions by ~5% on 
average, although is some cases the fraction is much larger (see 
Table|2]i. Due to the size of our original sample and the pseudo- 
random selection of H n-regions that are affected by these issues, 
we consider that this last cleaning process will have little effect 
in the overall statistical significance of our survey. 

Once derived the emission line intensities, we estimate their 
corresponding equivalent width for each H n region and line. For 



doing so, instead of using the classical procedure (i.e., measur- 
ing the flux within a narrow-band wavelength range centred in 
the line and in two adjacent ones corresponding to the contin- 
uum), we make use of the results from our fitting analysis. We 
derive the equivalent width by dividing the emission line inte- 
grated intensities by the underlysing continuum flux density. We 
estimated the continuum as the median intensity in a band-width 
of 100A, centred in the line, using the gas-subtracted spectra 
provided by our fitting procedure. With this method we can es- 
timate the equivalent width of nearby lines, which contaminate 
the measurements of this parameter using the classical method. 

Table illustrates the result of this analysis. It shows, for 
a sub-set of Hn regions in a particular galaxy (UGC 9837) the 
HB line intensity and relative flux of some of the most promi- 
nent emission lines. Table |6]reports the equivalent widths for the 
corresponding emission lines and regions. The same parameters 
are derived for all the H n regions in the different galaxies, as 
indicated in Appendix lAl 

6.1 .2. Structural parameters of the galaxies 

To understand the fundamental properties of the H n regions and 
their relation with the overall evolution of galaxies, it is re- 
quired to characterize the main structural parameters of these 
galaxies. We have collected the available information in public 
collections like NED (http : //ned. ipac . caltech.edu/) and 
Hyperleda (http://leda.univ-lyonl . fr/). Table Q] already 
contains the most relevant parameters for the current study, in- 
cluding the morphological type, the redshift, the integrated B- 
band magnitude and the B-V color. All the galaxies in the sam- 
ple are spirals by selection, but different kinds of spiral galaxies 
are covered, including galaxies with and without bars, galaxies 
showing rings, etc. The observed B-band magnitude is ~14 mag, 
in average. However, the galaxies selected from the PINGS sam- 
ple are in general brighter (B ~11 mag). As expected, galaxies 
have blue colors (B-V ~0.8 mag), although with a consider- 
able dispersion. The covered absolute magnitudes range from 

My 23 mag to 18 mag. In summary, the sample covers 

typical members of the so-called blue cloud, from typical stan- 
dard spiral galaxies to almost dwarfs. 
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Fig. 7. Left panel: [Ora] /15007/Hj8 vs. [Nn] /16583/Ha diagnostic diagram for the average properties of the Hn regions galaxy 
by galaxy listed in Table [7] The error bars indicate the standard deviation from the mean value. Middle panel: Similar diagnostic 
diagram for all the emission line regions detected by the described analysis with good quality measurements of both parameters 
(2230 regions). The color image and contours show the density distribution of these regions. The first contou r indicates the mean 
density, with each consecutive one increasing by four times this mean value. The [Kauffman n et al.l (120031) (red solid-line) and 
iKewlev et alj d2001l) (blue dashed-line) demarcation curves are usually invoked to distinguish between star-forming regions (below 
the red solid-line), and other source of ionization, like AGN/shocks/post-AGB (above the blue line). Regions between both lines 
are considered intermediate ones, indicating a mixed origin for the ionization. The error-bars at the top-left indicate the typical 
(mean) errors for the considered line ratios. Right panel: Similar diagram, including only the 124 regions at the core of each galaxy 
(r < 0.5r e , i.e., at ~ 0.25r f in average). The fraction of regions in the intermediate location is clearly higher. 



The listed information was complemented with additional 
parameters, like the maximum rotation velocity, the inclination, 
the position angle and the effective radius (defined as the radius 
at which one half of the total light of the system is emitted), 
derived from the analysis of the data presented here. The first 
three parameters were derived from the modeling of the gas ve- 
locity pattern extracted from the the Ha emission line fitting for 
the Hn regions, described in previous sections. The wide spatial 
coverage and high S/N of the Ha emission line in the integrated 
spectra for each region guarantee a good determination of the ve- 
locity pattern. The rotation curve w as fitted using a simple arctan 
model dStavelev-Smith et al.lll990h . 

2 Vi-ot 

v(r) = v sys + arctan(i • r - c) (3) 

n 

where v sxs is the systemic velocity of the gas and v rot is the 
asymptotic rotation speed of the disc, s characterizes the slope 
of v(r) in the inner part of the galaxy, r is the distance to the 
rotational center, and c is the parameter that characterizes any 
offset in the rotation axis of the galaxy. A model of the velocity 
map was created by re-projecting the best fitting arctan func- 
tion, taking into account the position angle and inclination of the 
galaxy. We fitted to the data following a similar procedure as the 
one described in Sa nchez et al.l (l2012 | ). using a ^-m inimization 
algorithm included in FIT3D dSanchez et al.ll2006bl) . 

As an initial guess for the fitting, the position angle and incli- 
nation were derived from the isophotal analysis described later 
in this section, and the maximum rotational velocity was set to 
half of the maximum difference in velocity from receding to ap- 
proaching velocities. For all the galaxies the rotational center 
is fixed to the location of the peak intensity in the V-band im- 
age created from the datacubes. The parameter c is fixed to zero 
(i.e., it is assumed that there is no offset between the rotation and 
photometric centers). The v sys is fixed to the median value of the 
gas velocities for those H n regions located in the inner regions 
(r <0.5r max , where r max is the maximum distance to the center 
for all the H n regions). 



Finally v 10t and s are fitted, together with the position angle 
and inclination of the galaxies, v ro , is fitted within a range be- 
tween ~0.3 and 6 times the maximum velocity difference among 
the Hn regions, and s is fitted between 0.1 and 10 arcsec 1 . Fig. 
[6]illustrates a typical result of this analysis, showing, for a partic- 
ular galaxy (UGC 9837) the Ha velocity map, the best fit model, 
and the residual. Despite the low inclination of the galaxies, in 
most of the cases it is possible to obtain a good model. In most 
of the cases the residual velocities ranges between ±15 km s , 
~15% of the maximum rotational velocity. This is expected due 
to random motions in the galax ies, compare e.g. lAndersen et ail 
( 20081); iNeumaver et all (1201 ll) . 

The effective radius was derived based on an analysis of 
the azimuthal surface brightness (SB) profile, derived based on 
elliptical isophotal fitting of the ancillary g-band images col- 
lecte d for the galaxie s (extracted from the SDSS imaging sur- 
vey, |%rkeUiU|200fl and Paper I). When these ancillary images 
were not available we used the B-band (Paper I). In order to 
homogenize the dataset, both sets of SBs were transformed to 
the B-band using the average g - B color for each galaxy. When 
both band images where available a comparison between the di- 
rectly derived and the estimated surface brightness profile was 
performed, finding no significant differences in the average gra- 
dient. We note here that the observed B and g-bands sample a 
range of wavelengths between ~A 1 50-4750A and -4450-52 1 OA, 
respectively, due to the redshift range of the sample. Therefore, 
there is an inherent imprecision in the intrinsic wavelength range 
in this analysis. 

The surface brightness profile was then fitted using a pure 
exponential profile, following the classical formula, 

/ = / exp[-(r/r d )] (4) 

where Ip is the central intensity, and rj is the disk scale-length 
(lFreemanlll970l) . using a simple polynomial regression fitting. 
Prior to this analysis, a visual inspection is performed to remove 
the inner-most values of the SB profile, strongly affected by see- 
ing effect, and/or not following a linear relation due to the pres- 
ence of other components like the bulge and/or bars. 
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The scale-length is used to derive the effective radius, de- 
fined as the radius at which the integrated flux is half of the total 
one, by integrating the previous formula, and deriving the rela- 
tion: 

r e = 1.67835r d (5) 

The results of these analyses are included in Table Q] In average, 
the derived inclination agrees with the visual selection of the 
galaxies as face-on spirals. The average inclination is ~33°, and 
only two galaxies have an inclination larger than 60° (NGC 7570 
and UGC5100). This confirms our visual classification as face- 
on galaxies. The average maximum rotational velocity is ~100 
km s _1 , with a wide range of values, between ~50 km s _1 and 
~300 km s" 1 (val ues which are typical for spiral galaxies, e.g., 
iPersic et al . 1996). Note that the effective radius ranges between 
~1.5 and ~5.5 kpc. 

7. Analysis and results 

In this section we analyze both the mean statistical properties of 
the H ii regions and explore the possible regular patterns in their 
radial variations. 

7.1. Statistical properties of the H n regions 

Despite the many different spectroscopic studies in extragalac- 
tic H ii regions, we still do not have the understanding of which 
are the statistical spectroscopic properties of these common star- 
forming regions. This is a fundamental problem that it is mostly 
due to the lack of big statistical samples, and the reduced number 
of coherent compilations. The lack of a well defined set of nor- 
mal values for the most frequent parameters, like the diagnostic 
line ratios (e.g., [Om]/Hy8 and/or [Nn]/Ha), ionization strength, 
dust attenuation and/or electron density is a clear limitation to 
understand if a particular set of H n regions is different from the 
average, and if different at which significance level. To address 
this question a statistically significant, large sample of Hn re- 
gions is required, with well derived spectroscopic parameters, 
over a large sample of star-forming galaxies of different types. 
An additional requirement is good spatial coverage, not biased 
towards the outer (bright) H n regions, which is a common bias 
in this kind of studies. Despite the large number of Hn regions 
catalogued in this work, the current sample is still incomplete to 
address this fundamental question. We will require a sample as 
the one that will be provided by a survey like CALIFA (Sanchez 
et al. 2012), without selection effects by galaxy types. However, 
the current catalogue of H n regions is good enough to derive the 
statistical properties of these regions for a sub-set of galaxies: 
quiescent/non highly disturbed, field, average luminosity spiral 
galaxies. 

Among the several different spectroscopic parameters de- 
scribing H ii regions, we have analyzed a set of them based on 
the strongest detected emission lines: (i) EWhq,, the equivalent 
width of the Ha emission line. This parameter is directly re- 
lated to the fraction of very young stars (~ 10 Myr), and can be 
used to estimate the aging process of the ionizing population; (ii) 
Ay, the dust attenuation derived from the Ha/H/3 Balmer decre- 
ment. To derive it the extinction law bv ICardelli et al.l <[l989) 
was assumed, with Ry =3.1, and the theoretical value for the 
unobscured line ratio for case B recombin ation of Ha/Hff = 
2.86, for r e =10,000K and « e =100cnT 3 (lOsterbrockl fl989): 
(iii) two typical diagnostic line ratios, [Om] /15007/H/? and 
[Nn] /16583/Ho', that define the nature of the ionization source; 



(iv) the ionization parameter, estimated as \og w U = -3.02 - 

0. 8010gi n ([Oii1/[On i1), a measureme nt of the strength of the 
ionization radiation (iDfaz et al.l [2000); (v) the oxygen abun- 
da nce, 12 + log(Q/H), de rived using the 03N2 indicator defined 
bv lPettini & Pagell d2004). that estimates the g as enrichment; (vi ) 
rim, the Stromgren radius of the H n regions (Osterbrock 1989), 

1. e., a hint of the size of the H ii regions based on pressure equi- 
librium considerations, and (vii) n e , the electron density derived 
from the [S n] doubl et line ratio, i.e., a proxy to the density of 
the ionized gas (e.g., Osterbrock 1989). 

Table [7] shows the mean values and standard deviations of 
the considered parameters for the H n regions of each galaxy, to- 
gether with the mean value for all the galaxies. Based on this 
analysis, it is possible to describe two kind of behaviors: on one 
hand, some of the parameters have a large variation object by 
object (e.g., EW[Hq,], Thu), reflecting the different physical con- 
ditions of Hn regions in individual galaxies. On the other hand, 
there are parameters with a well defined mean value and little 
variation object-by-object, and even region by region (e.g., the 
dust attenuation, the ionization parameter and the oxygen abun- 
dance), i.e., despite of the many differences between the consid- 
ered galaxies (luminosity, morphology, color), and the physical 
conditions in each H n region, it is possible to define a statisti- 
cally meaningful standard mean value for certain spectroscopic 
parameters. 

Due to the particular sample of galaxies studied here, and 
the large number of H n regions explored, we consider that these 
values define the average physical conditions of H ii regions for 
spiral galaxies in the Local Universe. They can either be used 
to determine whether a particular (spiral) galaxy deviates from 
the average population (i.e., it is metal rich or poor and/or it is 
more or less dusty), or as the anchor point of chemical evolu- 
tion of ionized gas along cosmological times. The Hn regions 
discussed here have a range of Ho- luminosities between 10 37 
and 10 41 erg s _1 , with a well defined bell-like shape centred at 
jq39.5±o.6 er g s -i Therefore, most of our regions correspond to 
intermediate/luminous ones, as indicated in Section[2] These H ii 
regions are expected to have a typical size of few to several hun- 
dreds of parsecs, at the edge of our spatial resolution or below. 

For some parameters the mean value is well defined (i.e., it 
shows a dispersion around the mean value of the order of the 
estimated error). Thus it is a good characterization of the con- 
sidered property of the ionized gas. However, it is important to 
remember that the variations within the considered distributions 
reflect changes in the physical conditions of the ionized nebu- 
lae, either within each galaxy, or galaxy by galaxy. To illustrate 
this effect we present in Fig. [7] thr ee panels showing the classi- 
cal [Om]/Hy6 vs. [Nn]/Ha BPT (iBaldwin etal.|[T9 81) diagnostic 
diagram for (i) the average values shown in Table [7j (ii) all the 
emission line regions detected using the procedure described in 
previous sections; and (iii) those emission line regions located at 
the center of the galaxies (r < Q.5r e ). In each panel, we include 
the Kauffmann et al. (2003) and Kewley et al. (2001) demar- 
cation curves. These curves are usually invoked to distinguish 
between star-forming regions, and other sources of ionization, 
like AGN/Shocks/post-AGBs. The location within both curves is 
normally assigned to a mix origin for the ionization and/or con- 
tamination by different sources of ionization. Since our sample 
is dominated by H n regions, most of the regions lie in the de- 
marcation region corresponding to star-forming areas, although 
a few of them are located in the so-called intermediate region. 
There are clear differences galaxy by galaxy, reflecting changes 
in the average oxygen abundance and ionization strength. On 
the other hand, it is clear that for the regions located in the cen- 
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Table 7. Median physical parameters of the ionized gas derived for the considered galaxies 
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/in i 


± {J.J 


-34.3 


± 7.9 


1.0 


± 


n a 


-U.DU = 


- U.zo 


-U.4z 


± 


n no 
U.U5 


-3.00 d 


= U.Z4 


O. /o 


± U.U/ 


T O 
Z.U = 


1 1 
Z 1.1 






line 1 99^n 


jy.O 


± yj.j 


-21.0 


± 3.7 


i .— 




yj.j 


n zL1 - 

-U.H-l _ 


. n in 


-yj.j i 




n OS 






O. / J 


-t- n 

=t= U.UJ 










UGC 12391 


39.7 


± 0.3 


-21.5 


± 4.7 


1.4 


± 


0.4 


-0.37 = 


= 0.18 


-0.48 


± 


0.06 






8.70 


±0.07 






3.9 + 


21.8 


NGC 628 


38.7 


±0.5 


-50.4 + 


26.7 


1.1 


± 


0.5 


-0.40 = 


= 0.28 


-0.55 


± 


0.09 


-3.75 d 


= 0.29 


8.68 


±0.10 


0.5 = 


= 0.5 


0.5 


±3.1 


NGC 1058 


38.7 


± 0.5 


-24.8 + 


11.7 


1.0 


± 


0.4 


-0.33 = 


= 0.33 


-0.50 


± 


0.12 


-3.82 d 


= 0.27 


8.70 


±0.12 


0.5 : 


= 0.3 


0.4 


± 2.3 


NGC 1637 


38.8 


± 0.5 


-19.0 


± 6.7 


1.3 


± 


0.4 


-0.52 = 


= 0.29 


-0.39 


± 


0.10 


-3.84 d 


= 0.25 


8.79 


±0.08 


0.5 : 


= 0.8 


1.5 


± 2.6 


NGC 3 184 


39.2 


±0.4 


-43.0 ± 


15.7 


1.3 


± 


0.5 


-0.64 = 


= 0.29 


-0.49 


± 


0.07 


-3.79 d 


= 0.27 


8.80 


±0.08 


1.0 = 


= 0.6 


0.7 


± 2.7 


NGC 3310 


39.8 


±0.8 


-133.0 + 


51.8 


0.2 


± 


0.2 


0.26 = 


= 0.12 


-0.70 


± 


0.12 


-3.31 d 


= 0.44 


8.42 


±0.07 


1.0 = 


= 1.4 


0.5 


± 1.8 


NGC 4625 


39.3 


±0.4 


-24.6 


±6.0 


0.6 


± 


0.4 


-0.55 = 


= 0.15 


-0.51 


± 


0.04 


-3.81 d 


= 0.16 


8.75 


±0.05 


1.0 = 


= 0.3 


0.3 


± 1.8 


NGC 5474 


38.4 


± 0.4 


-64.6 ± 


28.0 


0.7 


± 


0.4 


0.17 = 


= 0.20 


-0.91 


± 


0.10 


-3.49 d 


= 0.44 


8.39 


±0.08 


0.3 = 


= 0.2 


0.2 


± 3.0 


Mean 


39.5 


±0.6 


-28.4 + 


20.8 


1.3 


± 


0.4 


-0.25 = 


= 0.23 


-0.47 


± 


0.13 


-3.71 d 


= 0.18 


8.67 


± 0.10 


1.9 = 


= 1.8 


1.9 


±2.3 



Notes: The table shows the median values and standard deviations of the following physical parameters: (1) Decimal logarithm of the dust corrected 
luminosity of Ha in erg s _1 ; (2) Equivalent width of the pure Ho- emission line expressed in A(values decontaminated from the underlying stellar 
absorption); (3) Dust attenuation derived from the YLaftif} line ratio; (4) [OIII]/15007/H/3 line ratio, in decimal logarithm scale; (5) [NII]/16583/Hq' 
line ratio, in decimal logarithm scale; (6) Ionization parameter; (7) Oxygen abundance, 12 + log(0/H), derived using the 03N2 indicator; (8) 
Stomgren radius calculated for the HII regions in units of 100 parsecs; (9) Electron density in units of 100 electrons per cm 3 . 



ter of the galaxies the fraction of emission line regions lying in 
the intermediate region is larger. This is expected since these are 
the regions more likely to be contaminated by other ionization 
sources than star-formation: e.g., central low-intensity AGNs, 
Shocks, post-AGB stars. It is beyond the scope of this article to 
determine the real nature of this ionization in the central regions. 
This will require a more detail analysis and the use of other more 
suitable diagnostic diagrams like [NIIl/Ha vs. [SII]/Ha and/or 
[NIIl/Ha vs. [SII]6717,6731/i/a, to distinguish between HII re- 
gions, SNR and PNe. We just want to emphasize that if there is 
any possible contamination, this will more likely affect only the 
central regions. This clearly indicates a change in the ioniza- 
tion conditions galaxy by galaxy and across the optical extension 
of each galaxy. One may wonder whether the described differ- 
ences are induced by the errors in the represented parameters. 



However, most probably this is not the case. The typical error 
is represented as an error-bar in the central panel of the figure, 
showing that it is smaller than the dispersion of values found (in 
particular for the range of [Nn]/Ha values sampled). 

Another possible explanation is that the large range of red- 
shifts/distance has an effect in variation of some properties due 
to our coarse resolution. This could be the case for the differ- 
ences found for the Stomgren radius or the luminosity of Ha 
for some galaxies, like NGC 628, at a distance of ~9 Mpc and 
Mrk 1477, at a distance of ~90 Mpc. However, this cannot be the 
only explanation, since there are other galaxies with large differ- 
ence in both parameters at lower relative cosmological distance, 
like UGC 6410 (at -80 Mpc), and NGC 4109 (at -100 Mpc). 

In summary, although the mean values listed in Table [7] are 
somehow representative of the average ionization conditions in 
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log (R23) log (R23) log (03N2) 

Fig. 8. Left-panel: Distribution of the 03N2 vs. R23 line ratios for all the 1142 Hn regions with detected [On] A3121 emission 
line. The image and contours show the density distribution of both parameters. The first contours is at the mean density, with a 
regular spacing of four times this value for each consecutive contour. The red-dashed line shows the best fit found using a 2nd- 
order polynomial function between the two considered parameters. The light-blue solid-circles indicate the median values of both 
parameters, with their corresponding standard deviations represented as error bars, for consecutive bins of 0.10 in A03N2. The 
black solid li ne shows the relation expected between both parameter s when they are used to derive the same oxygen abundance, 
assuming the Trem onti et al.l (120041) fit for R23 and the lPetti ni & Pagel (2004) calibration for 03N2. The dot-dashed blue l ine s hows 
the same relation, if the oxygen abundance derived based on R23 was 70% of the one derived using the iTremonti et alJ d2004l) fit. 
Central panel: Similar plot for the distributions of the N202 line ratio vs. R23, for the same Hn regions. In this case the light-blue 
solid-circles indicate the median values and corresponding standard deviations of the parameters for consecutive bins of 0.10 in 
AN202. Right panel: Similar plot for the distribution of the N2 line ratio vs. 03N2. The blue solid circles indicate the median 
values of both parameters, with their corresponding standard deviations, for consecutive bins of 0.10 in A03N2. The red dashed- 
line shows the best fit found using a 2nd-order polynomial function between both parameters. In the three panels the typical/median 
errors of the represented parameters are represented as a black error-bar. 



each galaxy, they do not show the complete picture, since there 
are clear/expected variations of the ionizing condition across the 
optical extension of the galaxies. 



7.1 .1 . Strong-line calibrators of the oxygen abundance 

One of the fundamental parameters that can be derived from the 
spectroscopic analysis of Hn regions is the oxygen abundance. 
Oxygen is one of the easiest elements to measure in H n regions, 
due to the strength of its emission lines in the optical wavelength 
range. This is fortunate, since O is an a-process element made 
directly in short-lived massive stars (dominant in H 11 regions). It 
is a good proxy of all heavy elements, comprising ~50% of all 
the metals by mass in all the Universe. Therefore, it is a funda- 
mental parameter to understand the evolution of the stellar pop- 
ulations galaxy-by-galaxy and at different locations within the 
same galaxy. 

Accurate abundance measurements for the ionized gas 
in galaxies require the determination of the electron tem- 
perature (T e ) in this gas which is usually obtained from 
the ratio of auroral to nebula r line intensities , such as 
[Ora] ^4363/[Om] zl/14959,5007 dOsterbrocklll989l) . It is well 
known that this procedure is difficult to carry out for metal- 
rich galaxies since, as the metallicity increases, the electron 
temperature decreases (as the cooling is via metal lines) and 
the auroral lines eventually become too faint to measure. 
Therefore, calibrator s based on strong emission lines are used. 
lLopez-Sanchez et"a"fl d2012l) has recently presented a revision on 
the different methods, showing their main problems, and illus- 
trating in which range/conditions then can be applied. In partic- 
ular they describe the most widely used methods in large galaxy 
surveys, that rely on the measurement of different strong emis- 
sion lines (and line ratios) in the H 11 region spectrum and empir- 
ical calibrations with regions of well-known oxygen abundance. 



One of the first s trong emission line methods was proposed 
bv lPagel et all dl979t) . It relies upon the ratio of [O n] ,13727 and 
[Om] ,U4959,5007 with respect to Hy8, the so-called R23 ratio: 



R23 = 



7([Q 11] ,13727) + 7([Q ill] ,U4959, 5007) 



(6) 



However, this index is double valued, with two different calibra- 
tions for low metallicity (12 + log(0/H)<8.1) and high metal- 
licity (12 + log(0/H)>8.4) Hn regions. There is an ill-defined 
regime where regions with the same R23 (> 0.7) value have oxy- 
gen abundances that differ by almost and order of magnitude 
(e.g., see Figure Al in Lopez-Sanchez & Esteban 2010). Ratios 
such as Q3N2 dAlloin et al.lll979b iPettini & Pagell 2004]). N2Q2 



ftuu iii ci ai.iu?/?UJ~cuim ex. r atici 

(van Zee el al .l 119981: Dopita et al.l 120001), N2 dyan Zee et all 



119981: lPettini& Pagell 120041: iDenicolo et al l 120021) . and many 
others, were introduced in an attempt to solve this ambiguity in 
the derived abundances: 



03N2 
N202 
N2 



/([Om] A5007)/I(H/3) 
/([Nil] i6584)//(Ha) 
/([Nil] ,16584) 
/([On] ,13727) 
/([Nil] ,16584) 
I(Ha) 



(7) 
(8) 
(9) 



Some of these indicators are strongly affected by the dust at- 
tenuation, like N202 and R23, and should be determined af- 
ter correcting it. Others are less sensitive to dust attenuation. 
Furthermore, empirical calibrations based on direct estimates of 
the electron temperature of the ionized gas systematically pro- 
vide oxygen abundances which are systematically 0.2-0.4 dex 
lower than those de rived using calibrations bas ed on photoion- 
ization models (e.g. Lop ez-Sanche z et al. 20T2I). 
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To derive the main statistical properties of the oxygen abun- 
dance, we adopted the 03N2 ratio. This ratio is less depen- 
dent on the adopted correction for the dust attenuation, and 
uses emission lines covered by our wavelength range for all 
the galaxies in the current sample. Like the N202 indicator, it 
relies on the use of the intensity of [Nn]. In the high metal - 
licity regime -12+log(0/H)>8.3- nitrogen is produced by both 
massive and intermediate-mass stars (e.g. iPilvugin et all l2003h 
and hence nitrogen is essentially a secondary element in this 
metallicity regime. The [Nii]/Ha ratio will therefore become 
stronger with on-going star formation and evolution of galaxies, 
until very high metallicities are reached, i.e. 12 + log(O/H)>9.0, 
where [N n] starts to become weaker due to the very low electron 
temperature caused from strong cooling by metal ions. However, 
N2 is an indirect estimator that may depend on the ionization 
strength. 

We should remember here that for ~5Q% of our targets 
[O n] /13727 was not covered by the observed wavelength range, 
and therefore, the fraction of H n regions with R23 and N202 
line ratios is just half of the total number of regions, since both 
rely on this emission line. The direct estimate based on the elec- 
tron temperature is of less use. The fraction of H n regions with 
reliable measurements of the the auroral [Oin] /14363 is just 
~10%. This line is either too faint or is blended and/or strongly 
affected by the HgI4348 atmospheric light pollution emission 
line, and even Hy, in most of the cases. Due to that we con- 
sider 03N2 the best estimator for the purpose of the current 
study. In general it is considered that this estim ator is valid fo r 
the range between -1< log03N2 <1.9 (e.g. lYin et al]|2007l) . 
However it becomes les s reliable for m e tallicities lower than 
12 + log(0/H)<8.3 (e.g. lYin et al.l 120071: iLonez-Sanchez et alJ 
l20T2h . 

It is beyond the scope of the current study to make a 
detailed comparison of the oxygen abundances derived using 
the different p roposed methods , in a si milar way as it was 
presen ted by iKewlev & Ellison] d2008l) : lLopez-Sanchez et ail 
(120121) . However, it is a good sanity check to compare between 
some of the most frequently used ones. This comparison will 
help us to understand (1) if the H n regions described here are in 
the valid range for the adopted metallicity indicator; (2) which 
are the real differences in using one indicator or another; and (3) 
how accurate is the adopted correction for the dust attenuation, 
since we can compare indicators with different degree of depen- 
dency on the adopted correction. 

Fig. [8] shows the comparison between 03N2, N202 and N2 
indicators, all of them using the [Nn] /16584 line, with respect 
to the solely oxygen-based R23 ratio. Only the H n regions with 
detected [0 11] /13727 emission line have been included in this 
plot (1 142 individual regions). In each panel it is shown the im- 
age and contours of the density distribution of both parameters. 
The red-dashed lines show the best fits found using a 2nd-order 
polynomial function between the two considered parameters. 
The light-blue solid-circles indicate the median values of both 
parameters, with their corresponding standard deviations repre- 
sented as error bars, for consecutive equally spaced bins of 0. 10 
dex of the parameters shown. In all cases the distributions are 
similar to the ones previously described in the literature between 
the considered estimators an d/or th eir derived abundan ce (e.g . 
Ivan Zee et al.ll 19981: lYin et"aill2007l;lRo"sales-Ortega et alj201 lb . 
However, in none of the precedent studies the distributions are 
shown with such a large statistical sample. 

In the three cases there is a clear correlation between the me- 
dian/high dense distribution of both parameters. The correlation 
is stronger for N202 vs. R23 and for N2 vs. 03N2 than for the 




log (Mass lum /M ) 

Fig. 9. Relation between the stellar mass, in units of the so- 
lar masses, and the oxygen abundance at the central regions 
(r < 0.2 r e ), using the 03N2 calibration, for the galaxies con- 
sidered in the current study ( solid circles). The bla ck-solid line 
shows the relation found by Tre monti et alJ d2004l) for the star- 
forming galaxies in the SDSS sample, corrected for the mis- 
match between their suggested calibrator and the currently used 
here (i.e., -0.15 dex). The grey-dashed lines show the ~lcr and 
~2cr range around this relation. 



remaining one. The first of this relations is well described by a 
2rd order polynomial function: 

logN202 = 0.15 - 0.97 logR 23 - 0.39[logR 23 ] 2 (10) 

The dispersions of each parameter around the best fit curve 
are crmoi ~0.11 and cr R2 3 ~0.17, respectively. The dispersion 
in the y-axis is very similar to the typical (mean) errors of the 
corresponding line ratio, eN202 ~ 0.07. This indicate that most 
of the dispersion is due to the errors in the derived line ratios, 
and not in a physical dependence with a third parameter. 

The second relation is also well described by a 2rd order 
polynomical function: 

logN2 = 0.46-0. 12 log 03N2 - 0.19[logO3N2] 2 (11) 

For this distribution, the dispersions of each parameter 
around the best fit curve are cr w2 ~0.10 and ct 3n 2 ~0.32, re- 
spectively. Once more, the dispersion in the y-axis is very sim- 
ilar to the typical (mean) errors of the corresponding line ratio, 
<?N2 ~ 0.17. Again, this indicates that most of the dispersion is 
due to the errors in the measured line ratios, rather than a un- 
known relation with a third parameter. 

Finally, the correlation is less tight for the the 03N2 and R 2 3 
parameters, with a broader dispersion and a tail for R 2 3 > 0.7. In 
this case a 2nd order polynomial function describe will enough 

log03N2 = 0.18 -3.76 logR 23 +5.88[logR 23 ] 2 (12) 

However, the dispersion around the best curve is larger for the 
03N2 parameter, <xo3N2 ~0.62, but similar for the R 2 3 one, 
o"R23 ~0.17. This dispersion cannot be explained only by the 
larger typical (mean) error of the represented line ratios (eo3N2 ~ 
0.2 and eR 2 3 ~ 0.3, respectively), and most probably indicates a 
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real difference and/or a dependence with a third parameter (like 
the ionization strength). 

iTremonti et al.l ([2004) proposed a non linear fit of the oxygen 
abundance for the R23 indicator that is nowadays widely used. 
Assuming that this calibration derives a similar oxy gen abun- 
dance than the one proposed bv lPettini & Pagell d2004l) is is pos- 
sible to derive the following relation between both parameters: 



y = -1.406 + 0.987x + 0.825x 2 + 1.396x 3 



(13) 



where: y — log03N2 and x = logR23. This relation has 
a similar shape than the one described before for the range 
of values corresponding to 0.5 < R23 <1, however it is off- 
set by ~0.15 dex in R23 and ~ -0.5 dex in 03N2. The oxy- 
gen abundance d erived using both in dicators agrees if it is 
assumed that the ITremonti et alj d2004l) relation overestimate 
the abundance by ~30%. The differ ence is due to the pho- 



toionizati on models adopted by Tremonti et al.l d2004l) . as dis- 
cussed in Lopez-Sanch ez et al.l d2012l) . Correcting for that dif- 



ference, the derived relation between 03N2 and R23 agrees 
with the one found with our dataset (blue dashed-dotted line 
in Fig. [p. Therefore, the oxygen abundance der ived using the 
iPettini & Page! d2004l) and lTremonti et alJ d2004l) calibrators do 
not provide consistent results, as alread y pointed out in previous 
studies (e.g. lLopez-Sanchez et al.ll2012l) . However, the discrep- 
ancy is in the zero-point of the adopted calibration, rather than in 
the actual values for the indicators. The dispersion is tighter for 
R23 < 0.7, i.e., the parameter range for which the oxygen abun- 
dance is not double valued. For values of R23 larger than 0.7 
we found that the discrepancy with the derived relation shows a 
trend with the ionization parameter, being larger at lower values 
of XL 

iKewlev & Ellison! d2008l) performed a comparison between 
the oxygen abundance derived using different strong line indica- 
tors, based on analysis of large sample of SDSS spectra of star- 
forming galaxies. They derived a set of transformation for the 
different analyzed indicators, including R23 and 03N2. These 
transformation are consistent with the ones presented here, for 
the range of abundance for which their are valid (see Table 3 in 
IKewlev & Ellisonll2008l) . 

In summary, it seems that (i) N202 and R23 indicators pro- 
vide statistically similar empirical information on the oxygen 
abundance, i.e. measuring one or the other in the considered 
range of parameters will yield a similar estimation of the oxy- 
gen abundance, when a consistent calibration is used. On the 
other hand, 03N2 provides different information than the other 
two indicators, when dealing with individual regions, even when 
consistent calibrators are used; (ii) the range of values covered 
by the 03N2 indicator for the H 11 regions of our sample is within 
the values where this indicator is reliable (as indicated before); 
and (iii) the lack of outliers in the N202 vs. R23 distribution 
indicates that we have adopted a good correction for the dust 
attenuation. 



7.1.2. Mass-Metallicity Relation 

Stellar mass and metallicity are two of the most fundamental 
physical properties of galaxies. Both are metrics of the galaxy 
evolution process, the former reflecting the amount of gas locked 
up into stars, and the later reflecting the gas reprocessed by stars 
and any exchange of gas between the galaxy and its environ- 
ment. Understanding how these quantities evolve with time and 
in relation to one another is central to understanding the physical 
processes that govern the efficiency and timing of star formation 
in galaxies. 



ITremonti et alJ d2004l) found a tight correlation between stel- 
lar mass and gas-phase metallicity (MZR) spanning over 3 or- 
ders of magnitude in stellar mass and a factor of 10 in metallic- 
ity. This correlation has been interpreted as a direct evidence of 
the anticorrelation between metal loss with baryonic mass, and 
the ubiquity of galactic winds and their efficiency in removing 
metals from low-mass galaxies. Additional interpretations con- 
sider that the observed mass-metallicity relation is due to a more 
general relatio n between stellar, gas-phase metallicity and s tar 
formation rate dLara-Lopez et alJl2010tlMannucci et alJl20ldh . 

Despite the strength of the considered correlation and how 
it matches with our current understanding on galaxy evolution, 
there are a few caveats o n how the MZR is deri ved: (i) The oxy- 
gen abundance used by ITremonti et al.l d2004l) was derived as- 
suming a Bayesian approach based on simultaneous fits to all the 
most prominent lines assuming the results provided by photoion- 
ization models. However, as explained before these methods sys- 
tematically overestimate by 0.2-0.4 dex th e oxygen abundances 
derive d using the Te-based calibrations dLopez-Sanchez et al.l 
|2012|) ; (ii) the calibrator was derived for H n regions, but applied 
to integrated apertures in galaxies. Although it is widely used, it 
is is not firmly tested if this procedure is valid; (iii) the apertures 
of the considered spectra cover a different optical extent for each 
galaxy, due to the fixed aperture of the SDSS fiber spectra (3")- 

So far, there has been no major effort to test the reported rela- 
tion using spatially resolved spectroscopic information, like the 
one provided by the current analysis. The number/statistics of 
our analysis, galaxy-by-galaxy, is reduced, and we can only test 
if our results are consistent with the derived relation. However, it 
will be a first step towards a stronger statistical result, that could 
be provided with the la rge sample currently under observation 
by the CALIFA survey dSanchez et alj|2011l) . 

We derive the luminosity mass of the considered galaxies us- 
ing the integrated B-band magnitudes, B - V co lors and the av- 
erage mass-luminosity ratio (M/L) described in Bell & de Jond 
(2001 [). This derivation is less accurate than the one described 
in Tremonti et al . (2004). However, it should be enough for the 
purpose of this analysis, consid ering that the derived M/L has an 
expected scatter of ~0.15 dex dBell & de Jondl2001l) . As an ad- 
ditional cross-check we also derived the dynamical mass. We use 
the effective radius (r e ), the maximum rotational velocity (v rot ) 
listed in TableQ] and the classical formulaQ: 



Ma- 



ya 



(14) 



We note here that this mass is indeed the dynamical mass at a 
particular radius (roughly twice the effective one), and therefore, 
the derived values should be taken with caution for any further 
interpretation. As expected the dynamical mass is slightly larger 
than the stellar mass, since the first considers the total mass of 
the galaxy. However, both masses agree within +0.2 dex for 
~90% of the objects. Similar results have been recently reported 
for the same sort of galax ies, using more soph isticated deriva- 
tions of the stellar masses dGerssen et al. 2012). This consistent 
result strengthens the validity of our derived masses. 

Figure [9] shows the distribution of the oxygen abundance at 
the center of the galaxies vs. the stellar masses described before. 
The oxygen abundance has been derived extrapolating the radial 
gradients found for each galaxy (see Sec. I7.2K to the center. For 
comparison, we estimate the central abundance using the aver- 
age of the metallicity values within 0.2r e , using the 03N2 cal- 
ibrator, for those galaxies with H 11 regions at this distance (20 
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Table 8. Results of the analysis of the radial gradients of a set of properties. 



Parameter 


N gal 


'"cor 


o- r 


med a 


0"a 


err a 


med b 


en, 


err b 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 




25 


0.98 


0.02 


21.64 


0.84 


0.26 


1.68 


0.01 


0.23 


B-V 


25 


0.50 


0.22 


0.87 


0.20 


0.30 


-0.20 


0.15 


0.27 


ran (100 pc) 


11 


0.61 


0.24 


2.57 


0.25 


0.30 


-0.48 


0.20 


0.26 


12 + log 10 (OH) 


25 


0.60 


0.25 


8.84 


0.13 


0.30 


-0.12 


0.11 


0.27 


logwOEWaJ) 


25 


0.57 


0.21 


0.95 


0.32 


0.25 


0.51 


0.28 


0.22 


logi ([Nn]/Ha) 


25 


0.54 


0.24 


-0.41 


0.13 


0.30 


-0.10 


0.13 


0.27 


log 10 ([Oni]/H8 


25 


0.53 


0.28 


-0.59 


0.32 


0.25 


0.29 


0.26 


0.22 


log(f/) 


11 


0.35 


0.24 


-3.93 


0.29 


0.04 


0.09 


0.12 


0.08 


n E (100 cm 4 ) 


22 


0.21 


0.13 


2.48 


7.79 


0.26 


0.31 


0.67 


0.24 


A v (mag) 


25 


0.17 


0.17 


1.60 


0.88 


0.25 


-0.29 


0.52 


0.22 



Notes: (1) Parameter which gradient is analyzed, eg., PARAM = a + b * R/Re; (2) Number of galaxies fulfilling the criteria (n H n >30 regions); 
(3) median value of the derived correlation coefficients for the analyzed gradient; (4) standard deviation of the derived correlation coefficients; (5) 
median value for derived zero-points for each galaxy; (6) standard deviation of the derived zero-point; (7) median error of the derived zero-points; 
(8) median value of the derived slopes for each galaxy; (9) standard deviation of the derived slopes; (10) median error of the derived slopes. 



out of 38). Both values agree within the errors. We prefer to use 
the extrapolation of the derived gradient, since they have lower 
dispersion, and they can be derived for more galaxies (i.e., all 
but MRK 1477). We i nclude in Fig M the rel ation between both 
quantities derived by iTremonti et al. d2004l) . corrected for the 
mismatch between their calibrator and the currently used here 
(i.e., 0. 15 dex). The derived central abundance and stellar masses 
for the galaxies analyzed here are consistent with this relation, 
showing a similar distribution: about 75% of the galaxies are lo- 
cated at less than ~lcr of the scaled relation. 



processes is a possible flattening of the radial gradient and/or a 
truncation at outer regions, where some of these properties (e.g. 
oxygen abundance) show values that are expected at more in- 
ner regions. This change in the slope of the radial gradient of 
oxygen abund ance has been repo rted by several authors ( e.g. 
iBresolin et all 2009: Ma rino et al.ll2012l;lBresolin et alj|2012l) . If 
true, migration is an efficient mechanism to intermix the gas and 
stellar population pro perties of galaxies, and even to enrich the 
intergalactic medium dSanchez-Blazquez et alj|2009l for a dif- 
ferent explanation to the surface-brightness truncation). 



7.2. Radial gradients 

It is well known that different spectroscopic properties of Hn 
regions show strong variations across the area of disk galax- 
ies. In particular, some of these parameters (e.g., oxygen abun- 
dance, EW[Ha], etc.), show a strong radial gradient, that in 
average indicates that more evolved, metal rich, stellar popu- 
lations are located in the center of galaxie s, and less evolved , 
metal poor ones are i n the outer ones (e.g. IZaritskvet al.lll994t 
iRosales-O rtega et ai] |201 U and references therein). This obser- 
vational result is consistent with the our current understand- 
ing of the formation and evolution of spiral galaxies (e.g. 
Tsuiimot o et al.ll2~010i , and references therein). Gas accretion 
brings gas into the inner region, where it first reaches the re- 
quired density to ignite star formation. Thus the inner regions 
are populated by older stars, and they have suffered a faster gas 
repro cessing, and galaxies experience an inside-out mass growt h 
(e.g. iMatteucci & Francois! fl989t iBoissier & Prantzosl Tl999h . 
Both the extinction-correcte d color gradients in nearby galaxies, 
M unoz-Mateos et al.l (120071) . and weak dependence of the mass- 
size relation with redsh ift dTruiillo et al. 2004-1 iBarden et al.l 
120051; iTruiillo et alj|200ol) support an inside-out scenario for the 
evolution of disks. 

Pure inside-out growth is affected by internal and external 
events, like bars and/or galaxy interaction and mergers, which 
induce radial migration. In principle, radial migration affects the 
stellar populations more strongly, since they are dynamically 
more susceptible to the effects that produce migration (interac- 
tions, bars, etc.). The overall result is a change in the slope of 
the surface bright ness profile and colors of these galaxies (e.g. 
Bakos et al. 2008h. However, gas clumps , like Hn regions, may 
also be affected (e.g. lMinchev et alJ2012l) . The net effect of these 



Despite of the several different studies describing these ob- 
servational events, there is a large degree of discrepancy be- 
tween the actual derived parameters describing the gradients: 
(i) slope of the gradient; (ii) average value and dispersion of 
the zero-point and (iii) scale-length of the truncation. In gen- 
eral, this is mostly due to different observational biases and his- 
torical methods to perform the analysis: (i) most of the stud- 
ies lack a proper statistical number of analyzed H n regions per 
galaxy; (ii) in many cases the samples cover mostly the outer re- 
gions (R>R e ), with a lack of Hn regions in the inner ones; (iii) 
there is no uniform method to analyze the gradients. In some 
cases the phys ical scales of the g alaxies (i.e., the radii in kpc) 
are used, (e.g. lMarino et al.ll2.Q12h . In others the scale-length are 
normalized to the R25 radii, i.e., the radii at which the surface 
brightness in the ff-band re ach the value of 25 mag/arcsec 2 (e.g. 
IRosales-Ortega et"aT]l201 ll) . Finally, a much more reduced num- 
ber of studies tries to normalize the sc ale-length ba sed on the 
effective radii (described in Sec. 16.1.21 . Diaz (1989), in one of 
the earliest compilations of H n regions already showed that the 
effective radius is the best one to normalize the abundance gra- 
dients. This is somehow expected since it is directly related with 
the mass concentration and reflects how fast the gas is recycled 
in different galaxies. Despite its benefits, it is the least used one, 
which has produced an ill-defined set of definitions for these gra- 
dients. It is important to remember here that both the physical 
scale of the radial distance and/or the normalized one to an abso- 
lute parameter like the R25 radius, although they are widely used, 
could not produce useful gradients to compare galaxy to galaxy, 
since in both cases the derived gradient is correlated with either 
the scale-length of the galaxy of its absolute intensity. 
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Fig. 10. Summary of the analysis of the gradients of a set of physical properties of the H n regions in the considered galaxies. 
Each panel shows, from left to right, (i) the distribution of correlations coefficients of a each of the considered parameter (from 
top to bottom: the surface brightness, the oxygen abundance and the dust attenuation) with respect to the radial distance; (ii) the 
distribution of the zero-points of the linear regression and (iii) the distribution of the slope of the same regression. The orange solid 
line represents, for each of the last two histograms, the expected histogram in case of a Gaussian distribution of the data, assuming 
the mean and standard-deviation of the distribution of each analyzed parameter, and sampled with same bins. 



7.3. Analysis of gradients galaxy by galaxy 

We used our catalogue of H n regions to characterize the radial 
gradients of the physical properties of the Hn regions. For each 
galaxy we derive the correlation coefficient, and the slope and 
zero-point of a linear regression; for the radial distribution we 
derived each of the parameters discussed in the previous sec- 
tion. The radial distance is normalized to the effective radius of 
each galaxy, listed in Table [4] Only those galaxies with a good 



sampling of the radial distribution of H n regions have been con- 
sidered. To do so, we adopted two criteria: (i) the galaxy should 
comprises more than 30 individual Hn regions, and (ii) there 
should be equally distributed along the sampled region (~2.5 

Re). 

Table[8]summarizes the result of this analysis. It includes, for 
each of the considered parameters, the actual number of galax- 
ies considered(A^ Sfl /), the median of the correlation coefficients 
(r cor ), with its standard deviation (o>), the median zero-point 
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Fig. 11. Left panel: Radial distribution of the equivalent width of Ha (in logarithm scale of its absolute value), once scaled to the 
average value at the effective radius for each galaxy. The image and contours show the density distribution of H n regions in this 
parameter space. The first contour indicates the mean density, with a regular spacing of four times this value for each consecutive 
contour. The light-blue solid-circles indicate the mean value (with its corresponding standard deviation represented as error bars), 
for each consecutive radial bin of ~0.15 R e . The solid line shows the average linear regression found for each individual galaxy, 
as listed in Table [8] The red dashed line shows the actual regression found for all the Hn regions detected for all the galaxies. The 
average error of the equivalent width of Ha for the H n regions represented in this plot is shown by a single error bar located at 
the top-left side of the panel. Right panel: Similar radial distribution for the oxygen abundance derived using the 03N2 indicator, 
once scaled to the average value at the effective radius for each galaxy, following the same format. The average error of the derived 
oxygen abundance (without considering systematic errors) is shown by a single error bar located a the top-right side of the panel. 
The right-side scale shows the corresponding value for the 03N2 indicator. The solid-orange square indicate the average abundance 
of the solar neighborhood, at the distance of the Sun to the Milky-Way galactic center (in effective radius). 



(med a ), with its standard deviation (cr a ), and the mean error of 
each individual zero-point (err a ), together with the median slope 
(medb), its standard deviation (a-/,), and the mean error of each 
individual derived slope (err/,). In addition to the physical pa- 
rameters of the H n regions discussed in the previous section we 
have included, for comparison purposes two properties of the un- 
derlying stellar population, that are well known to correlate with 
the radial distance: the Z?-band surface brightness (/UbX an d the 
B — V color. The table is ordered by the strength of the correlation 
for the physical parameters of the H n regions. 

The overall analysis is illustrated is Fig. [TO] where we show 
the distribution of the correlation coefficients, zero-points, and 
slopes for each individual galaxy, for three particular parame- 
ters: (i) the surface brightness, as an example of a strongly cor- 
related parameter; (ii) the oxygen abundance, i.e., the strongest 
correlated parameter for the ionized gas; and (iii) the dust ex- 
tinction, a parameter with no gradient in most of the cases. For 
highly correlated parameters the mean value, shown in Table [7] 
is not enough to characterize their behavior. 

There are parameters that show a statistically significant cor- 
relation with the radial distance (r cor >0.5, and <j cor <0.3). These 
parameters are the Stromgren radii, the oxygen abundance, the 
equivalent width of Ho- and the classical ionization diagnos- 
tic line rations ([Oiii]/H/? and [Nn]/Har). The first correlation 
may indicate that the H n regions are larger towards the center. 
However, we have to be cautious about this statement, since we 
do not have a direct estimate of the actual size of the H n regions, 
and the Stromgre n radii are based o n particular radiation equilib- 
rium conditions (Osterbrock 1989). We need to recall here that 
since many of our so-called H n regions, are indeed, H n aggre- 
gates, this radius should be considered as an effective radius. The 
next two correlations show that the more metal rich and evolved 
stellar populations are located in the inner regions, in agreement 



with expectations. The remaining two gradients indicates that in 
general the strength of the ionization is larger in the outer than 
in the inner regions. Finally, the ionization parameter shows a 
weaker trend, in the way that less powerful H n regions are lo- 
cated in the central areas, while the stronger ones are located in 
the outer regions. However, this result has to be taken with care, 
since the correlations are weak and the ionization parameter de- 
pends both on the electron density and the number of ionizing 
photons, which may present also radial gradients. 

Finally, there are parameters with little dependence with the 
radial distance, like the dust attenuation and the electron density. 
If the radial gradient of the parameters discussed before (e.g., 
oxygen abundance and/or equivalent width of Ha) is a conse- 
quence of the time evolution of the H n regions (aging, enrich- 
ment, lose of strength and expansion), those later ones are not 
strongly affected by this overall evolution. 

7.4. Universal gradients 

For those properties showing a strong correlation, it is equally 
important to understand if the gradient is universal, within our 
range of explored parameters. Inspecting the histograms shown 
in Fig.fTOland the values for the slopes listed in Table[8] we can 
conclude that all the parameters showing either a clear correla- 
tion or a trend show well defined slopes, taking into account the 
considered errors for the measured parameters. For the equiv- 
alent width of Ho- and the oxygen abundance a Lilliefors-test 
indicates that the histogram of the slop es of the gradien ts are 
consistent with a Gaussian distribution (Lilliefors 1967). This 
implies that we can define a characteristic value for the slope 
and that we do not find a population of galaxies with slopes in- 
consistent with the normal distribution. To illustrate this result 
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we show in Fig. [TT] the radial distribution of two of the param- 
eters with stronger correlation: (i) the equivalent width of Ha 
and (ii) the oxygen abundance. We include all the H n regions 
in our catalogue (~2600). For doing so, a global offset has been 
applied galaxy by galaxy to normalize the equivalent width and 
the abundance at effective radius. This offset was derived by sub- 
tracting the value derived by the individual regression to the me- 
dian value listed in Table|8] A new regression is estimated for all 
the dataset of Hn regions. Both the average regression listed in 
Table |8] (red dashed-line), and the regression for the full sample 
(black solid-line), has been included. It is clear that both cor- 
relations are very similar, which support our claim that there is 
a universal slope for the gradient of either the equivalent width 
and the oxygen abundance when the radial scale is normalized 
by the effective radius. 

Fig. [TT] shows that the scatter around the described corre- 
lation, ~0.6 dex, is much larger than the average error of each 
individual measurement, ~0.15 dex, for the case of the equiva- 
lent width of Ha. However, in the case of the oxygen abundances 
both values are very similar, ~0.2 dex. If we interpret this scatter 
as a typical coherent scale-length, we can conclude that stellar 
populations with similar star-formation histories are distributed 
at a wider radial ranges than the corresponding ionized gas with 
similar oxygen abundance. Thus, the possible intermix due to ra- 
dial migrations (or similar effects) affects strongly the stars than 
the ionized gas, as expected. A larger relative scatter in the prop- 
erties of the stellar population compared with th e oxygen a bun- 
dance ones was reported even in our Galaxy (e.g. lFriell ll995). 

For the particular case of the oxygen abundance, the same 
caveats expressed in Sect. 17.1.11 have to be taken into account 
her e. We derive the abund ance using the 03N2 empirical estima- 
tor dPettini & Pagelll2004 . calibration which has a certain range 
of validity. Due to that we included the corresponding scale of 
03N2 values to illustrate that we are in the range where this cal- 
ibrator provides reliable results. We repeated the analysis using 
the R23 calibrator, deriving similar results for the radial gradi- 
ent. However, there is a larger scatter in this distribution, mostly 
due to the larger error bar of our derivation of oxygen abun- 
dance based on this later method, and the lower number of H n 
regions with measured [On] A3127 emission line (~50% of the 
total sample). 

Despite these caveats, the main result are valid: i.e., there 
seems to be a universal radial gradient for oxygen abundance 
when normalized with the effective radii of the galaxies. This re- 
sult apparently contradicts previous ones, in which the slope of 
the radial gradients show a trend on certain morphologica l char- 
acteristics of the galaxies. lVira-Costas & Edmunds! dl992h shows 
that barred spirals have a shallower gradient than non-barred 
ones. However, we must recall here that this statements are 
based, in general, in gradients on physical scales (i.e. dex kpc -1 ), 
not in normalized ones. A slight decre ase in the scatter between 
the differen t slop es is appreciated by IVila-Costas & Edmunds! 
(fl99l and |5ia3 0989), when the scale-length is normalized 
by the disk-scale, which agrees with our results. We conjec- 
ture that both results would come into agreement if the size- 
luminosity/morphological type relation was considered. 

In order to make an independent cross-check we analyzed 
the distribution of gradients among our galaxies for barred (10 
galaxies) and un-barred ones (28 galaxies). No significant dif- 
ference is found between the distributions, or between any of 
them and the one comprising the full sample of galaxies; on the 
basis of a f-test analysis, the probability of being different was 
63.1% in the worst case. In addition to this test we repeated 
the comparison for those galaxies with and without well de- 



fined spiral arms (see Table [3}. The nature of flocculent spirals 
is still unclear, although the star formation in these galaxies (and 
therefore, the H 11 regions) may not be connected with density 
waves, but rather wi th differentially rotating galactic disks (e.g. 
lEfremov et al.l fl989). They could be produced also due to gravi- 
tation al instabilities in the stars due to interactions (e.g. lToomrd 
11964 . We found no differences among the distributions of slopes 
for both kind of galaxies. 

Finally, for the Grand Design galaxies (i.e., those galaxies 
with two clearly defined spiral arms listed in Table 0), we re- 
peated the analysis of the gradients for the Hn regions located 
at each spiral arm. Despite of the lower number statistics, we 
found not significant difference between the gradients for each 
spiral arm. 

Fig. QT| shows that for the spatial scale which is well sam- 
pled by our catalogue of Hn regions (r< 2r<,), the mean value 
of the analyzed properties follow the described distribution. For 
radii larger than ~2 effective radii, the mean abundance and Ha 
equivalent width seem to be slightly offset, with values compat- 
ible with more inner regions (~ 1 r e ). Although this seem to 
agree with previous results reporting a flattening on the oxygen 
abundance, and even with the migration as a possible explana- 
tion, we have to take this result with a lot of caution. There is 
just a hand-full of H 11 regions sampled by our catalogue in this 
spatial regime (less than a 1%), and therefore the result is not 
significant. On the other hand, it seems that there is a decrease 
of mean oxygen abundance in the inner regions of the galaxies 
(<0.3 r e ), compared with the expected value extrapolating the 
reported correlation. This decrease is not due to the abundance 
indicator used, since we detect it using either 03N2 or R23 one. 
Indeed, this decrease was already detected in individual galax- 
ies, li ke in NGC 628 dSanchez et al.ll201 UlRosales-Ortega et all 

Eoub. 

We explored the possibility of a statistical origin of this ef- 
fect, by analyzing the distribution of oxygen abundance galaxy 
by galaxy. However, we found that this decrease is present in 20 
of the 28 galaxies with H 11 regions detected in the central region. 
A lower slope in the abundance gradient may induce this effect. 
However, we found no significant difference in their gradients, 
with respect to the average values. 

Finally, we have explored which would be the actual location 
of the solar neighborhood in the proposed abundance gradient 
shown in Fig.Qj] For the distance of the Sun t o the center of the 
Milky way we adopted the value deduced by Eisenhauer et alj 
(120031) . of 7.94+0.42 kpc. The value of the disk scale-length, 
and therefore the effective radius, is still controversial, with large 
discrepancies between different results. The more likely value is 
~4.5 kpc (i.e., r e ~ 7.6 kpc), although a dispersion larger than 
1 kpc is expected (Ivan der Kruit & Freeman! 1201 ll) . Therefore 
the Sun would be at ~1.1 effective radius of the Galactic cen- 
ter. We note th at the abundanc e in the solar neighborhood is also 
under debate dAsplu nd 2007). The current conve ntion is that it 
is aro und ~8.7, with a dispersion of ~0.2 dex ( Asplun d et alj 
2006). Adopting this scale-length and abundance the actual lo- 
cation of the solar neighborhood would be just at the average 
location from our derive radial gradient (orange square in Fig. 
rXTT t. Considering that the oxygen abundance has been scaled to 
the average value at the effective radius, it is interesting to note 
that our galaxy behaves as the average for typical spiral galaxies. 

8. Summary and Conclusions 

In this article we have analyzed the H 11 regions of a sample of 
38 nearby galaxies. The sample has been constructed selecting 
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the almost face-on spiral galaxies presented in Paper I, and sim- 
ilar on es extracted from the PINGS survev lRosales-Ortega et al] 
d2010t) . In both cases, the galaxies have been observed using 
IFS, using the same instrument and sampling similar wavelength 
ranges with similar resolutions. The IFS data cover most of the 
optical extension of these galaxies, up to ~2-2.5 effective radii, 
with enough S/N for the purpose of this study. 

We presented HIIexlorer, an automatic method to detect and 
segregate Hn regions/aggregations based on the contrast of the 
Ha intensity maps extracted from IFS datacubes. Once detected, 
the procedure provides us with a segmentation map that is used 
to extract the integrated spectra of each individual H n region. 
The method has been compared successfully with other codes 
available in the literature for NGC628, the largest galaxy (in 
projected size) included in the sample. 

In total, we have detected 3 107 H n regions, 2573 all of them 
with good spectroscopic information (i.e. ~60 Hn regions per 
galaxy). This is by far the largest nearby, 2-dimensional spectro- 
scopic survey presented for this kind of regions up-to-date. Even 
more, our selection criteria and dataset guarantee that we cover 
the galaxies in an unbiased way, regarding the spatial sampling. 
Therefore, the final sample is well suited to understand the spa- 
tial variation of the spectroscopic properties of H n regions in 
this kind of galaxies. 

A well-tested automatic decoupling procedure has been ap- 
plied to remove the underlying stellar population, and to de- 
rive the main properties (intensity, dispersion and velocity) of 
the strongest emission lines in the considering wavelength range 
(covering from [On] /13727 to [S n] /16731). A final catalogue of 
the spectroscopic properties of these regions has been created for 
each galaxy. Additional information regarding the morphology, 
spiral structure, gas kinematics, and surface brightness of the un- 
derlying stellar population has been added to each catalogue. In 
particular, a detailed analysis of the number and structure of the 
spiral arms has been presented, associating each H n region to 
the nearest spiral arm, and classifying them among inter-arms 
and intra-arms ones. This analysis is highly experimental, and 
its main aim is to understand asymmetries in the spiral distribu- 
tion of H ii regions. To our knowledge this is the first time that an 
H ii region catalogue is provided with two dimensional informa- 
tion and linked to the spiral structure and intra-arm association. 
The full capabilities of the derived catalogues will be presented 
in subsequent articles. 

In the current study we focused on the understanding of the 
average properties of the H n regions and their radial distribu- 
tions. Among the analyzed properties we include: (1) the equiv- 
alent width of Ha, (2) the amount of dust attenuation (Ay), (3) 
the most widely used diagnostic line ratios (^p and -^r), (4) 
the ionizing parameter (U), (5) the oxygen abundance (O/H), (6) 
the effective Stromgren radius and (7) the electron density. 

We derive the following results: 

1 . There is statistically significant change in the ionization con- 
ditions across the field of view. The fraction of H n regions 
affected either by AGNs, shocks or post-AGB stars is much 
larger for the central regions (r < 0.25r e ), than in the outer 
ones. 

2. The average abundances for the central regions are com- 
parable with t he mass-metallicity relation derived by 
iTremontie t al. (2004]). This is particularly important con- 
sidering that both the technique to derive the stellar masses 
and the central metallicities (including the adopted calibra- 
tor) are different. The fully independent confirmation of the 
proposed relation through integral field spectroscopy is par- 



ticularly important as aperture biases can be virtually elimi- 
nated. 

3. The oxygen abundance and the equivalent width of Ha 
present a radial gradient that, statistically, has the same slope 
for all the galaxies in our sample, when normalized to the 
effective radius. The derived slopes for each galaxy are com- 
patible with a Gaussian random distribution and are indepen- 
dent of the morphology of the analyzed galaxies (barred/non- 
barred, grand-design/flocculent). 

4. Other properties show no gradient but seem to vary across 
each galaxy, and galaxy by galaxy (like the electron den- 
sity), without a clear characteristic value, or they are well 
described by the average value either galaxy by galaxy or 
among the different galaxies (like the dust attenuation). 

The Hn regions catalogues of this analysis will 
be made publicly available for a better use on 
the astronomical community (see Appendix [A] and 
|f tp : //ftp . caha . es/CA LIFA/e arly_studies/HI I/tables/ 
)• 

Acknowledgements. We thank the referee for the detailed and helpful comments 
that have improved considerably the quality of this article. 

We thank the director of CEFCA, Dr. M. Moles, for his sincere support to 
this project. 

We thank the Viabilidad , Diseno , Acceso y Mejora funding program, 
ICTS-2009-10, for funding the data acquisition of this project. 

S.F.S., F.F.R.O. and D. Mast thank the Plan National de Investigation y 
Desarrollo funding programs, AYA2010-221 1 1-C03-03 and AYA2010-10904E, 
of the Spanish Ministerio de Ciencia e Innovation, for the support given to this 
project. S.F.S. wants to note that the project continues despite of the lack of 
support of the same program during the 2012 period. 

F.F.R.O. acknowledges the Mexican National Council for Science and 
Technology (CONACYT) for financial support under the programme Estancias 
Posdoctorales y SabAaticas al Extranjero para la ConsolidaciA§n de Grupos de 
InvestigaciAsn, 2010-2011 

J.M. and J. P. acknowledge financial support from the Spanish grant 
AYA2010-15169 and Junta de Andalucia TIC114 and Excellence Project P08- 
TIC-03531. 

D. M. and A. M.-I. are supported by the Spanish Research Council within 
the program JAE-Doc, Junta para la Ampliation de Estudios, co-funded by the 
FSE. 

RA. Marino was also funded by the Spanish programme of International 
Campus of Excellence Moncloa (CEI). 

J. I. -P., J. M. V., A. M.-I. and C. K. have been partially funded by the projects 
AYA2010-21887 from the Spanish PNAYA, CSD2006 - 00070 "1st Science with 
GTC" from the CONSOLIDER 2010 programme of the Spanish MICINN, and 
TIC1 14 Galaxias y Cosmologfa of the Junta de Andalucia (Spain). 

This paper makes use of the Sloan Digital Sky Survey data. Funding for the 
SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the 
Participating Institutions, the National Science Foundation, the U.S. Department 
of Energy, the National Aeronautics and Space Administration, the Japanese 
Monbukagakusho, the Max Planck Society, and the Higher Education Funding 
Council for England. The SDSS Web Site is http://www.sdss.org/ 

The SDSS is managed by the Astrophysical Research Consortium for 
the Participating Institutions. The Participating Institutions are the American 
Museum of Natural History, Astrophysical Institute Potsdam, University of 
Basel, University of Cambridge, Case Western Reserve University, University 
of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, 
the Japan Participation Group, Johns Hopkins University, the Joint Institute 
for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and 
Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences 
(LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for 
Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New 
Mexico State University, Ohio State University, University of Pittsburgh, 
University of Portsmouth, Princeton University, the United States Naval 
Observatory, and the University of Washington. 

References 

Alloin, D., Collin-Souffrin, S., Joly, M., & Vigroux, L. 1979, A&A, 78, 200 
Andersen, D. R., Walcher, C. J., Boker, T., et al. 2008, ApJ, 688, 990 



24 



S. F. Sanchez et al.: Integral field spectroscopy of a sample of nearby galaxies: 



Asplund, M. 2007, in IAU Symposium, Vol. 239, IAU Symposium, ed. 

T. Kuroda, H. Sugama, R. Kanno, & M. Okamoto, 122-129 
Asplund, M., Grevesse, N., & Sauval, A. J. 2006, Communications in 

Asteroseismology, 147, 76 
Bakos, J., Trujillo, I., & Pohlen, M. 2008, ApJ, 683, L103 
Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 
Barden, M., Rix, H.-W., Somerville, R. S., et al. 2005, ApJ, 635, 959 
Bell, E. F. & de Jong, R. S. 2001, ApJ, 550, 212 
Bertin, E. & Amouts, S. 1996, A&AS, 117, 393 
Boissier, S. & Prantzos, N. 1999, MNRAS, 307, 857 

Bradley, T. R., Knapen, J. H., Beckman, J. E., & Folkes, S. L. 2006, A&A, 459, 
L13 

Bresolin, E, Kennicutt, R. C, & Ryan-Weber, E. 2012, ArXiv e-prints 
Bresolin, E, Ryan- Weber, E.. Kennicutt, R. C, & Goddard, Q. 2009, ApJ, 695, 
580 

Cappellari, M. & Emsellem, E. 2004, PASP, 116, 138 

Cardelli, J. A., Clayton, G. C, & Mathis, J. S. 1989, ApJ, 345, 245 

Cid Femandes, R., Mateus, A., Sodre, L., Stasiriska, G, & Gomes, J. M. 2005, 

MNRAS, 358, 363 
Corwin, Jr., H. G, Buta, R. J., & de Vaucouleurs, G. 1994, AJ, 108, 2128 
Denicolo, G, Terlevich, R., & Terlevich, E. 2002, MNRAS, 330, 69 
Diaz, A. I. 1989, in Evolutionary Phenomena in Galaxies, ed. J. E. Beckman & 

B. E. J. Pagel, 377-397 
Diaz, A. I., Castellanos, M., Terlevich, E., & Luisa Garcia- Vargas, M. 2000, 

MNRAS, 318, 462 

Dopita, M. A., Kewley, L. J., Heisler, C. A., & Sutherland, R. S. 2000, ApJ, 542, 
224 

Efremov, Y. N., Korchagin, V. I., Marochnik, L. S., & Suchkov, A. A. 1989, 

Uspekhi Fizicheskikh Nauk, 157, 599 
Eisenhauer, E, Schodel, R., Genzel, R., et al. 2003, ApJ, 597, L121 
Elmegreen, D. M. & Elmegreen, B. G. 1987, ApJ, 314, 3 
Fathi, K., Beckman, J. E., Zurita, A., et al. 2007, A&A, 466, 905 
Freeman, K. C. 1970, ApJ, 160, 811 
Friel, E. D. 1995, ARA&A, 33, 381 

Garcia-Benito, R., Perez, E., Diaz, A. I., Maiz Apellaniz, J., & Cervino, M. 201 1, 

AJ, 141, 126 
Gamett, D. R. 2002, ApJ, 581, 1019 

Gavazzi, G, Zaccardo, A., Sanvito, G, Boselli, A., & Bonfanti, C. 2004, A&A, 
417, 499 

Gerssen, J., Wilman, D. J., & Christensen, L. 2012, MNRAS, 420, 197 
Gonzalez Delgado, R. M. & Perez, E. 1997, ApJS, 108, 199 
James, P. A., Shane, N. S., Beckman, J. E., et al. 2004, A&A, 414, 23 
James, P. A., Shane, N. S., Knapen, J. H., Etherton, J., & Percival, S. M. 2005, 
A&A, 429, 851 

Jansen, R. A. 2000, PhD thesis, Kapteyn Astronomical Institute, Groningen, The 
Netherlands, and Harvard-Smithsonian Center for Astrophysics, Cambridge, 
MA 

Kauffmann, G, Heckman, T. M., Tremonti, C, et al. 2003, MNRAS, 346, 1055 
Kehrig, C, Monreal-Ibero, A., Papaderos, P., et al. 2012, A&A, 540, All 
Kelz, A., Verheijen, M. A. W., Roth, M. M., et al. 2006, PASP, 118, 129 
Kennicutt, Jr., R. C. 1983, ApJ, 272, 54 
Kennicutt, Jr., R. C. 1992, ApJ, 388, 310 

Kennicutt, Jr., R. C, Armus, L., Bendo, G, et al. 2003, PASP, 115, 928 
Kennicutt, Jr., R. C. & Kent, S. M. 1983, AJ, 88, 1094 

Kewley, L. J., Dopita, M. A., Sutherland, R. S., Heisler, C. A., & Trevena, J. 

2001, ApJ, 556, 121 
Kewley, L. J. & Ellison, S. L. 2008, ApJ, 681, 1 183 
Koleva, M., Prugniel, P., Bouchard, A., & Wu, Y. 2009, A&A, 501, 1269 
Lara-Lopez, M. A., Bongiovanni, A., Cepa, J., et al. 2010, A&A, 519, A31 
Laurikainen, E., Salo, H., Buta, R., Knapen, J. H., & Comeron, S. 2010, 

MNRAS, 405, 1089 
Lequeux, J., Peimbert, M., Rayo, J. F., Serrano, A., & Torres-Peimbert, S. 1979, 

A&A, 80, 155 
Lilliefors, H. W. 1967, JASA, 62, 399 

Lopez, L. A., Krumholz, M. R., Bolatto, A. D., Prochaska, J. X., & Ramirez- 
Ruiz, E. 2011, ApJ, 731,91 
Lopez-Sanchez, A. R. 2010, A&A, 521, A63 

Lopez-Sanchez, A. R., Dopita, M. A., Kewley, L. J., et al. 2012, ArXiv e-prints 

Lopez-Sanchez, A. R. & Esteban, C. 2008, A&A, 491, 131 

Lopez-Sanchez, A. R. & Esteban, C. 2010, A&A, 517, A85 

MacArthur, L. A., Gonzalez, J. J., & Courteau, S. 2009, MNRAS, 395, 28 

Mannucci, E, Cresci, G, Maiolino, R., Marconi, A., & Gnerucci, A. 2010, 

MNRAS, 408, 2115 
Marino, R. A., Gil de Paz, A., Castillo-Morales, A., et al. 2012, ApJ, 754, 61 
Marmol-Queralto, E., Sanchez, S. E, Marino, R. A., et al. 201 1, A&A, 534, A8 
Martin, P. & Roy, J.-R. 1994, ApJ, 424, 599 
Matteucci, F. & Francois, P. 1989, MNRAS, 239, 885 
Minchev, I., Famaey, B., Quillen, A. C, et al. 2012, ArXiv e-prints 



Moustakas, J. & Kennicutt, Jr., R. C. 2006, ApJS, 164, 81 
Munoz-Mateos, J. C, Gil de Paz, A., Boissier, S., et al. 2007, ApJ, 658, 1006 
Neumayer, N, Walcher, C. J., Andersen, D., et al. 201 1, MNRAS, 413, 1875 
Ocvirk, P., Pichon, C, Lancon, A., & Thiebaut, E. 2006, MNRAS, 365, 46 
Oey, M. S., Parker, J. S., Mikles, V. J., & Zhang, X. 2003, AJ, 126, 2317 
Osterbrock, D. E. 1989, Astrophysics of gaseous nebulae and active galactic 

nuclei (University Science Books) 
Pagel, B. E. J., Edmunds, M. G, Blackwell, D. E., Chun, M. S., & Smith, G. 

1979, MNRAS, 189, 95 
Paturel, G, Petit, C, Prugniel, P., et al. 2003, A&A, 412, 45 
Perez-Gonzalez, P. G, Zamorano, J., Gallego, J., Aragon-Salamanca, A., & Gil 

de Paz, A. 2003, ApJ, 591, 827 
Persic, M., Salucci, P., & Stel, F. 1996, MNRAS, 281, 27 
Pettini, M. & Pagel, B. E. J. 2004, MNRAS, 348, L59 
Pilyugin, L. S., Thuan, T. X., & Vilchez, J. M. 2003, A&A, 397, 487 
Ringermacher, H. I. & Mead, L. R. 2009, MNRAS, 397, 164 
Rosales-Ortega, F. E, Diaz, A. I., Kennicutt, R. C, & Sanchez, S. F. 2011, 

MNRAS, 415, 2439 
Rosales-Ortega, F. E, Kennicutt, R. C, Sanchez, S. E, et al. 2010, MNRAS, 405, 

735 

Roth, M. M., Kelz, A., Fechner, T, et al. 2005, PASP, 1 17, 620 
Sanchez, S. F. 2006a, Astronomische Nachrichten, 327, 850 
Sanchez, S. E, Aceituno, J., Thiele, U., Perez-Ramirez, D., & Alves, J. 2007a, 
PASP, 119, 1186 

Sanchez, S. E, Cardiel, N, Verheijen, M. A. W., Pedraz, S., & Covone, G. 2007b, 
MNRAS, 376, 125 

Sanchez, S. E, Garcia-Lorenzo, B., Jahnke, K, et al. 2006b, New A Rev., 49, 
501 

Sanchez, S. E, Jahnke, K, Wisotzki, L., et al. 2004, ApJ, 614, 586 
Sanchez, S. E, Kennicutt, R. C, Gil de Paz, A., et al. 2012, A&A, 538, A8 
Sanchez, S. E, Rosales-Ortega, F. E, Kennicutt, R. C, et al. 201 1, MNRAS, 410, 
313 

Sanchez-Blazquez, P., Courty, S., Gibson, B. K, & Brook, C. B. 2009, MNRAS, 
398, 591 

Sarzi, M., Falcon-Barroso, J., Davies, R. L., et al. 2006, MNRAS, 366, 1151 
Skillman, E. D. 1989, ApJ, 347, 883 

Staveley-Smith, L., Bland, J., Axon, D. J., Davies, R. D., & Sharpies, R. M. 

1990, ApJ, 364, 23 
Thilker, D. A., Braun, R., & Walterbos, R. A. M. 2000, AJ, 120, 3070 
Toomre, A. 1964, ApJ, 139, 1217 

Tremonti, C. A., Heckman, T. M., Kauffmann, G, et al. 2004, ApJ, 613, 898 

Trujillo, I., Forster Schreiber, N. M., Rudnick, G, et al. 2006, ApJ, 650, 18 

Trujillo, I., Rudnick, G, Rix, H.-W., et al. 2004, ApJ, 604, 521 

Tsujimoto, T., Bland-Hawthorn, J., & Freeman, K. C. 2010, PASJ, 62, 447 

van der Kruit, P. C. & Freeman, K. C. 201 1, ARA&A, 49, 301 

van Zee, L., Salzer, J. J., Haynes, M. P., O'Donoghue, A. A., & Balonek, T. J. 

1998, AJ, 116, 2805 
Vazdekis, A., Sanchez-Blazquez, P., Falcon-Barroso, J., et al. 2010, MNRAS, 

404, 1639 

Viironen, K, Sanchez, S. E, Marmol-Queralto, E., et al. 2012, A&A, 538, A144 
Vila-Costas, M. B. & Edmunds, M. G. 1992, MNRAS, 259, 121 
Walcher, J., Groves, B., Budavari, T, & Dale, D. 2011, Ap&SS, 331, 1 
Yin, S. Y, Liang, Y. C, Hammer, E, et al. 2007, A&A, 462, 535 
York, D. G, Adelman, J., Anderson, Jr., J. E., et al. 2000, AJ, 120, 1579 
Zaritsky, D., Kennicutt, Jr., R. C, & Huchra, J. P. 1994, ApJ, 420, 87 



Appendix A: Catalogues of the H n regions 

The results of the overall analysis on the properties of the H n re- 
gions is compiled in three different catalogues per galaxy, each 
of them comprising different information. The nomenclature 
of each catalogue is table . HII . TYPE . GALNAME . csv (for the 
coma separated version) or table . HII . TYPE . GALNAME . txt 
(for the space separated ones), where TYPE correspond to each 
of the following types: coords (coordinates of the Hn regions), 
flux (fluxes of the stronger emission lines) or EW (equivalent 
width of the stronger emission lines). GALNAME corresponds to 
the galaxy name as listed in Table Q] 

The catalogues are stored in the CALIFA public ftp server: 
|f tp : //ftp . caha . es/CALIFA/early_studies/HIl7ta bles/ 
The content of each table is described in the following sections. 
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Fig. A.l. Left-panel: Spatial distribution of the [O m]/H/3 line ratio along the spatial extension of the UGC 108 1 galaxy, derived from 
the analysis of the Hn regions. Each symbol corresponds to a Hn region, its filling color corresponding to the shown parameter, 
scaled as displayed in the right-size color-table. Red symbols correspond to Hn regions associated with arm 1, and black ones 
correspond to those associated with arm 2 (where the indexing of the arms was selected in arbitrary way). Grey symbols represent 
H ii regions without a clear association with a particular arm, following the criteria described in Sect. [5] The circles represent those 
H ii regions below the iKauffmann etaD J2003) demarcation line in the BPT diagram shown in F ig. [71 and the squares corresponds 
to those ones located in the intermediate region between that line and the Kewl ev et alJ {2001) one. The size of the symbols are 
proportional to the Ha intensity. Right-panel: Similar spatial distribution for the absolute value of the equivalent width of Ha in 
logarithm scale. Symbols are similar to those described for the left panel. 



A. 1. coords tables 

They include all the information regarding the location of each 
H ii region within the galaxy, and additional information regard- 
ing their relation with the galaxy morphology, kinematics and 
the Ha luminosity. The tables have the following columns: 

1 . ID, unique identifier of the H n region, including the name 
of the galaxy (NAME) and a running index (NN), in the form: 
NAME-NNN. 

2. RA, the right ascension of the H n region. 

3. DEC, the declination of the H n region. 

4. Xobs, the relative distance in right ascension to the center of 
the galaxy, in arcsec. 

5. Yobs, the relative distance in declination to the center of the 
galaxy, in arcsec. 

6. Xres, the deprojected and derotated distance in the X-axis 
from the center, in kpc. 

7. Yres, the deprojected and derotated distance in the Y-axis 
from the center, in kpc. 

8. R, the deprojected and derotated distance to the center, in 
kpc. 

9. Theta, the deprojected position angle of the Hn region, in 
degrees. 

10. ^Arm, the ID of the nearest spiral arm. 

1 1 . Inter^ m , a flag indicating if the H n region is most probably 
associated to a particular arm, 1, or most likely an inter-arm 
region, 0. 

12. D m i n -Arm> the minimum distance in arcsec to the nearest spiral 
arm. 

13. T>Arm, the spiralcentric distance, i.e., the distance in arcsec 
along the nearest spiral arm. 



14. Theta^™, the angular distance in degree to the nearest spiral 
arm. 

15. vel/f ff , the Ha rotational velocity, in km s 

16. e_vel#Q., error of the Ho- rotational velocity, in km s _1 . 

17. loglOCLtfoO, decimal logarithm of the dust corrected lumi- 
nosity of Ha, in using of Erg s . 

18. e_logl0(L// ff ), error of the decimal logarithm of the dust 
corrected luminosity of Ha, in using of Erg s _1 . 

A.2. flux tables 

They comprise the fluxes and ratios of the stronger emission 
lines detected in each H ii region, as derived from the analysis 
described in Sect. 16. II They include the following columns: 

1 . ID, unique identifier of the H n region, described in the pre- 
vious section. 

2. FIuxh/3, observed flux of H/3, in units of 10~ 16 erg s _1 cm -2 . 

3. e_Flux///3, error of the observed flux of Hf3, in units of 10~ 16 
erg s _1 cirT 2 . 

4. rat_0II_Hb , [O n] A3721/H/3 line ratio. 

5. e_rat_0II_Hb, error of the [O n] XilllfHp line ratio. 

6. rat_0III_Hb, [Om] /15007/H/3 line ratio. 

7. e_rat_0III_Hb, error of the [Om] /15007/HyS line ratio er- 
ror. 

8. rat_0I_Hb, [O i] ^6300/Hj6 line ratio. 

9. e_rat_0I_Hb, error of the [O i] A63QQ/H/3 line ratio. 

10. rat_Ha_Hb, Ha/H/3 line ratio. 

1 1 . e_rat_Ha_Hb, error of the Ha/H/? line ratio. 

12. rat_NII_Hb, [Nil] /16583/H/3 line ratio. 

13. e_rat_NII_Hb, error of the [N n] /16583/H/3 line ratio. 

14. rat_He_Hb, Hel 6678/H/3 line ratio. 
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15. e_rat_He_Hb, error of the Hel 6678/H/3 line ratio. 

16. rat_SII 6717_Hb, [S n] ,167 17/H/J line ratio. 

17. e_rat_SII6717_Hb, error of the [Sn] ,167 17/H/3 line ratio. 

18. rat_SII6731_Hb, [Sn] M73 1/HJ3 line ratio. 

19. e_rat_SII6731_Hb, error of the [Sn] ,16731/H/J line ratio. 

20. BPT_type, flag indicating the location of the ionized gas re- 
gion in the classical BPT [Oiii]/H/3 vs. [Nn]/Ha diagnostic 
diagram shown in Fig. [7] with the following values: (0) loca- 
tion undetermined, due to the lack of any of the required line 
ratios with suffi cient S/N for this a nalysis; (1) region located 
below the Kauffman n et al. (2003) demarcation line, i.e., cor- 
responding to a classical star-forming H n region; (2) region 
locate d in the i ntermediate a r ea wit hin the lKauffmann et all 
d2003l) and the iKewlev et alj ( 20011) demarcation lines, i.e.; 
(3) region located above the IKewlev et alj (2001) demarca- 
tion line, i.e., corresponding to the area expected from being 
ionized by an AGN and/or a shock. 

A.3. EW tables 

They comprise the equivalent width, in Angstroms, of the 
stronger emission lines detected in each H n region, as derived 
from the analysis described in Sect. 16.11 They include the fol- 
lowing columns: 

1 . ID, unique identifier of the H n region, described in the pre- 
vious section. 

2. EW_0II, equivalent width of [On] ,13727. 

3. e_EW_0II, error of the equivalent width of [O n] ,13727. 

4. EW_Hbeta, equivalent width of the HB emission line. 

5. e_EW_Hbeta, error of the equivalent width of HB. 

6. EW_0III, equivalent width of [O m] ,15007. 

7. e_EW_0III, error of the equivalent width of [O in] ,15007. 

8. EW_0I, equivalent width of [O i] ,16300. 

9. e_EW_0I, error of the equivalent width of [O i] /16300. 

10. EW_Halpha, equivalent width of Ha. 

1 1 . e_EW_Halpha, error of the equivalent width of Ha. 

12. EW_NII, equivalent width of [Nn] ,16583. 

13. e_EW_NII, error of the equivalent width of [N n] ,16583. 

14. EW_SII, equivalent width of the [S n] ,1,16717,6731 doublet. 

15. e_EW_SII, error of the equivalent width of the 
[Sn] ,1,16717,6731 doublet. 

An example of the use of the parameters listed in these ta- 
bles is shown in Fig. [6] and in the analysis performed in Sect. 
16.1.21 Fig. IA.1I shows an example of the two dimensional dis- 
tribution of the properties included in the described catalogues, 
which makes use of the many different properties studied for the 
H ii regions. 

Appendix B: Empirical correction of the [N n] 
contamination in Ha narrow band images 

The Ha luminosity observed in spiral and irregular galaxies 
is believed to be a direct tracer of the ionization of the inter- 
stellar medium (ISM) by the ultraviolet (UV) radiation which is 
produced by young high-mass OB stars. Since only high-mass, 
short-lived, stars contribute significantly to the integrated ioniz- 
ing flux, this luminosity is a direct tracer of the star formation 
rate (SFR), independent of the previous star formation history. 
This is why dust-corrected Ha luminosity is one of the most 
widely used observables to derive SFR in galaxies. 

Among the several methods used to derive the Ha luminos- 
ity and its distribution across the optical extension of the galax- 



ies, narrow-band image is by far the cheapest one in terms of 
telescope time and complexity. Thus, it is the most frequentl y 
used too (e.g., I James et all 120041; iPerez-Gonzalez et al .1 12003). 
The technique is rather simple: (i) the galaxy is observed using 
both a narrow-band filter centered at the wavelength of Ha and 
a broader filter covering a wider range around the same wave- 
length range; (ii) The broad-band image is used to correct for the 
underlying continuum and provide a relative flux calibration of 
the continuum subtracted narrow-band image; (iii) by perform- 
ing a flux calibration of the broad-band image it is possible to 
have an accurate calibration of the decontaminated Ha emission 
map. 

One of the major limitations of narrow-band Ha imaging is 
contamination by the [Nn] line doublet, which ca n hardly be 
derive d from a single narrow-band imaging (e.g. Ijames et aU 
2005). Most frequently used narrow-band filters have a width 
of ~50-80A, and in most of the cases the emission from the 
[Nn] /U6548,6583 contaminate the derived Ha emission map, 
and therefore, the derived Ha luminosity and SFR. Despite of the 
fact that the [N n]/Ha line ratio present a strong variation a cross 
the field for star-forming galaxies (e.g. San chez et al.ll201ll and 
FigGK it is generally assumed an average correction. 

The most co mmonly used cor rect ions for entire galaxies ar e 
those derived bv lKennicutil (119831) and lKennicutt & Kentld 19831) . 
Spectrophotometric [N n]/Ha ratios of individual extragalactic 
Hn regions from the literature (see Kennicutt & Kent 1983 
and references therein) were compiled from 14 spiral galax- 
ies (mostly of type Sc) and 7 irregular galaxies. The average 
Ha/(Ha + [Nn]) ratio was found to be fairly constant, spanning 
the ranges 0.75±0.12 for the spirals, and 0.93+0.05 for the ir- 
regulars (Kennicutt 1983). In terms of the ratio [Nn] -total/Ha 
this corresponds to a median value of 0.33 for spirals and 0.08 
for irregulars. These values were calculated by finding the [N n]- 
total/Ha ratio of the brightest H n regions, averaging for each 
galaxy and then determining the mean value for spiral and irreg- 
ular types. This implicitly assumes that all Hn regions have the 
same proportion of [Nn]-total to Ha emission as those regions 
measured. We have already illustrated that this may not be the 
case in general, e.g., Fig. Q However, recent results show that 
the total integrated spectra of galaxies may have a stronger [N n] 
emission than the one reported before, with a larg e variation line 
rati os for differen t galaxy types Kennicutt (1992). 

iJansenl d2000l) already showed that there is a trend of the av- 
erage [N n]/Ha ratio with galaxy luminosity. In their study, the 
lines ratios for galaxies bright er than Mp = -1 9.5 are in agree- 
ment with the values found bv lKennicuttldl992l) . but fainter than 
this a striking trend is seen towa rds much lower valu es of this 
ratio. A similar trend is found bv lGavazzi et al.l ([2004). 

To our knowledge the only attempt to make a correction 
ac ross the optical ex tension of each galaxy is the one introduced 
bv lJames et al.l d2004l) . In this study they used the spatial distri- 
bution of the [N n]/Ha line ratio, instead of an average correction 
for all the entire galaxy. However, their estimations of the distri- 
bution of the line ratio are based on physical principles rather 
than in direct measurements. 

Making use of our extensive catalogue of individual emis- 
sion line regions, we have explored possible empirical correc- 
tions that allow us to perform spatially resolved corrections 
across the optical extension of each individual galaxy. For do- 
ing so we correlated the measured [N n]/Ha line ratios for each 
individual region with different parameters easy to address us- 
ing broad-band photometry, like the Ha intensity, B and V-band 
magnitudes, B - V colors and radial distance relative to the ef- 
fective radius. 
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Fig. B.l. left panel: Distribution of the [Nn]/Ha line ratio along the B-V color for all the detected Hn regions. The image and 
contours show the density distribution in this space of parameters. The first contour is at the mean density, with a regular spacing 
of four times this value for each consecutive contour. The solid line shows the actual regression found between the two represented 
parameters, right panel: Histograms of the relative difference between the Ha intensity derived from the narrow band images 
described in Sect. [3] and the ones derived from the fitting technique over the extracted spectra, described in Sect. 16.11 before (blue- 
solid histogram) and after (red-shaded histogram) applying the correction from the [N n] contamination based on the correlation 
shown in the left-panel. 



Among the different explored linear relations, the one with 
stronger correlation coefficient (r ~0.5), and better defined zero- 
point and slope, was the one with the B-V color: 

log 10 ([Nll]/#a) = -0.64 ±0 .2o + 0.36 ±0 .i 8 (fi - V) (B.l) 

Fig. IB. II left panel, shows the distribution of [Nn]/Ha line 
ratios along the B-V color for all the detected H n regions. The 
image and contours show the density distribution in this space 
of parameters, together with a solid line showing the best fitted 
linear regression described before. This relation can be used as a 
simple proxy of the [N n]/Ha line ratio, and used to decontami- 
nate the described Ha narrow-band images. 

We applied the proposed correction to the Ha narrow-band 
images described in Sect. [3] to demonstrate its improvement over 
a classical single correction over the entire galaxy. We derive 
for each detected H n region the corresponding flux in the Ha 
narrow-band images, using the segmentation maps provided by 
HIIexplorer (Ha lmg ). Then, we obtain the relative difference be- 
tween this contaminated flux and the real Ha flux measured us- 
ing the fitting procedure describe in Sect. 16.11 The mean value 
of this difference is ~0.28±0.19, which is consistent with the 
typical average correction found in the literature (e.g., Kennicutt 
1983). Fig. IB. II right panel, shows a solid-blue histogram of the 
relative difference between the Ha flux derived from the narrow- 
band images, and the ones derived using the emission line fitting 
procedure after applying this average correction. 

The derived correction is then applied on the same data using 
a iterative processes. In each iteration the B-V color is used to 
guess the [N n]/Ha line ratio. This ratio, together with the decon- 
taminated Ha flux derived from in the previous iteration, is used 
to derive the [N n] intensity. Finally, this intensity is subtracted to 
the original contaminated Ha lmg intensity to derive a new decon- 
taminated flux. For the first iteration it is assumed that the inten- 
sity decontaminated by the average correction is a good starting 
estimation of the Ha flux. After three iterations the decontami- 
nated Ha flux converge with a few percent. Fig. IB. II right panel, 
shows a hashed-red histogram of the relative difference between 



the the Ha flux derived from the narrow-band images, and the 
ones derived using the emission line fitting procedure after ap- 
plying the proposed correction. It is clear that the new histogram 
has a lower dispersion that the previous one, with a standard de- 
viation of ~0. 15. The proposed correction improves the accuracy 
of the derived Ha intensity by a a ~60%, in average. 

The effect of this correction is stronger when analyzing 
the spatial distribution of the Ha emission and/or the relative 
strength of the SFR. As we already indicated the [N n]/Ha ratio 
tends to decrease with the radius. Therefore, an average correc- 
tion would overestimate the Ha emission and the SFR in the in- 
ner regions, and underestimate it in the outer ones. Indeed, some 
authors have corrected the image-based Ha fluxes for [N n] con- 
tamination using available sp ectroscopical data of selected H n 
regions within the galaxy (e.g. lLopez-Sanchez & Esteban2008). 

So far, we did pay no attention to the possible physical con- 
nections between the explored parameters. The described rela- 
tion may indicate a physical connection between the conditions 
of the ionized gas and the ionizing stellar population, or between 
both of them and a third parameter not considered here. A pos- 
sible origin of this connection could be a co-evolution of both 
components of the Hn regions. The study of this connection, 
that we will address in future works, is beyond the scope of the 
current analysis. 
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